Algebraic integrable dynamical systems, 2+1-dimensional models in 
wholly discrete space-time, and inhomogeneous models in 
2-dimensional statistical physics 



I.G. Korepanov, 

Chelyabinsk University of technology, 
76 Lenin prosp., Chelyabinsk 454080, Russia 



Abstract 

This paper is devoted to constructing and studying exactly solvable dynamical systems in 
discrete time obtained from some algebraic operations on matrices, to reductions of such systems 
leading to classical field theory models in 2+1-dimensional wholly discrete space-time, and to 
connection between those field theories and inhomogoneous models in 2-dimensional statistical 
physics. 
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Chapter 1 



Dynamical system connected with 
the dimer model 



1. Definition of the dynamical system. Gauge invariance 

In this chapter we consider the following dynamical system in discrete time. Let 

L-( A B 

be a block matrix, A, . . . , D being n x n matrices consisting of complex numbers. Consider the 
following two operations: construction of the inverse matrix 



L=\ i I \^L*=< A C 



and the block transposing 

7 ~ ' C L) ) ~ " ~ I B L) 
Now let a (birational) mapping / be a composition of these two operations: 

f(L) = (L-y. (i.i) 

Let us introduce the discrete integer-valued time r, and let the matrix L depend on r so that 

L(r + l) = /(L(r)). (1.2) 

This dynamical system has been already mentioned in literature |3Qj . In the present chapter, 
the integrability of this system is demonstrated, assuming that the "motion" is considered up to 
a "gauge transformation" (see below). We also demonstrate what conditions (reduction) must be 
imposed on the matrix L in order to obtain a meaningful evolution equation in 2+1-dimensional 
space-time. This equation is a 2+1-dimensional version of Toda lattice in discrete time. The 
evolution is of hyperbolic character: perturbations propagate not faster than fixed "light speed" . 
The solution to a Cauchy problem can be constructed in theta functions according to a rather 
simple, in principle, scheme. The remarkable property of the 2+1-dimensional model is that its 
"integrals of motion" are nothing else than a statistical sum of the well-known flat dimer model 
(in our case, the statistical sum of dimer model depends on two "spectral parameters"). 



Returning then back to the general dynamical system (1.2), we consider its connections with a 



discrete analogue of the Lax L,A pair and possible generalizations coming therefrom. 

Definition 1.1. We will call the gauge transformation of the matrix L the following transforma- 
tion of its blocks: 

A -> GAff, ...,£)—> GDH, (1.3) 
with G and H non-degenerate n x n matrices. 
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It is clear that the transformation (TO) commutes with the shift by two time units. Two 



matrices L and L' connected by the transformation (1.3) will be called gauge equivalent. Thus, 



dynamics fll.2| ) induces a dynamics on the set of gauge invariance classes of matrices L. 

2. Vacuum curves and vacuum vectors 



It turns out that the dynamics ( |l.2| ) preserves the so-called vacuum curve F of the operator L (the 
bases being fixed, we make no difference between a linear operator and its matrix). To be exact, T 
remains unchanged under the transformation fof, and undergoes a simple transformation under /. 
The curve T together with the class of linear equivalence of the pole divisor of the vacuum vectors 
(see below) determines the matrix L up to a gauge transformation. The set of those classes of linear 



equivalence is isomorphic to a complex torus — the Jacobian of the curve T. The dynamics (1.2) 
linearizes on the Jacobian, i.e. the transformation / corresponds to a constant shift on the torus. 
Now, let us discuss these facts in detail. 

The vacuum curve of the operator L is an algebraic curve in the space C 2 of two variables u,v. 
Here are two equivalent definitions of it . 

Definition 2.1. Consider the relation 

L{U®X) = V®Y, (2.1) 

wherein 

"-(!)■ V ={1 

are two-dimensional vectors, X and Y are n-dimensional vectors. For a generic matrix L, the 



non-zero solutions (U, V, X, Y) of the relation (2.1) are parametrized, up to a scalar factor in X 
and Y, by points of an algebraic curve T of genus g = {n — l) 2 given by an equation of the form 

P(u,v) = 0, (2.2) 

P(u, v) being a polynomial of degree n in each variable, i.e. 

n 

P{u,v) = a jk u J v k . (2.3) 
F is called the vacuum curve of the operator L. 



Definition 2.2. The vacuum curve of the operator L is the curve T in C 2 given by the equation 
P(u, v) = detiV^LU) = det(uA + B — uvC - vD) = 0, (2.4) 

where 

V ± = (l,-v). 

Let us denote the points of the vacuum curve by the letter z = (u,v) € T. Then U = U{z) 
and V — V(z) are meromorphic vectors on T with the pole divisors Djj and Dy of degree n, while 
X = X(z) and Y = Y(z), if normalized by, e.g., the condition that their nth coordinates equal 
unity, become meromorphic vectors with pole divisors Dx and Dy of degree n 2 — n [[l2| . Under 
this normalization, a meromorphic scalar factor h(z) must be added into ( [2.l| ): 

L(U(z) ® X{z)) = h(z)V{z) ® Y{z). (2.5) 

The linear equivalence of divisors 

Du + Dx ~ Dy + Dy 
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holds and is provided by the function h(z) in the sense that h{z) has its poles in the points of 
Djj + Dx and zeros in the points of Dy + Dy. 

As is shown in the paper E^|, the vacuum curve equation P{u, v) = and the class of linear 
equivalence of divisor D x or Dy determine a generic matrix L to within a gauge transformation, 
and vice versa, the gauge transformations do not change the vacuum curve and the classes of linear 
equivalence of divisors. In other words, the correspondence 

(class of gauge equivalence of L) <-> (r, the class of D x ) 

is a birational isomorphism. 

We will call X{z) the vacuum vector and Y{z) the covacuum vector in the point z of the curve 
r. X(z) — X(u,v) generates the (one-dimensional) kernel of the matrix 

uA + B - uvC - vD. (2.6) 



The Definition 2.1 allows one to trace what happens with the vacuum curve and vacuum vectors 
under the transformation L — » L , while the Definition 2.2 allows one to trace what happens under 
the transformation L — > L . Namely, it is seen from the relation 

L~ l {V{z) <8> Y(z)) = h(z)- 1 U{z) ® X(z) 

that the vacuum curve equation for the matrix L _1 is 

P{v,u) = 0, 

while its vacuum vector in the point (v, u) coincides with the covacuum vector of the initial matrix 
L. As for the block transposing, the vacuum curve equation for the matrix L* 

det{uA + C - uvB - vD) = 

may be rewritten as 

u n v n detiv^A -B + u- x v- x C - u^D) = 0, 

The vacuum vector of the matrix L l in the point (— v , — u^ 1 ) of its vacuum curve coincides with 
the vacuum vector X(u, v) of the matrix L. 

Combining these considerations, one finds out that the vacuum curve T of the matrix (L^ 1 )* 
is given by equation 

u n v n P(~u-\-v- 1 ) = 0, 

while the vacuum vector X(— it , —v^ 1 ) coincides with the vector Y(u, v) of the matrix L. 
Identifying the curves T and T by means of the isomorphism 

(u,v) <-» (-it -1 , -v^ 1 ), 

one sees that 

D x ~D Y ~D X +Du-D v , 



i.e. 



which means that, in essence, the transformation (1.1) results in adding a fixed element of the 
Picard group, namely the equivalence class of the divisor Djj — Dy, to the pole divisor Dx of the 
vacuum vectors. It is clear also that after two transformations one returns to the initial curve: 

f = r. 
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Figure 1.1. Integrable dynamical system in the 2+1-dimensional space-time 

3. Reduction to evolution equation in the 2+1-dimensional 
space-time 

The dynamical system of the previous section admits an interesting reduction, i.e. some special 
choice of the matrices A, . . . , D that is in agreement with the evolution. In this section, it will be 
convenient to treat the matrices A, ...,£) as linear operators acting from the linear space Tii into 
the linear space Tii (of the same finite dimension). This being the situation at the moment r, the 
operators act, of course, from Ti.2 into Tit at the moment r + 1, and so on. 

Let each of the spaces Tit, Tii be a direct sum of lm/2 identical subspaces of dimension d, 
where I, m are even numbers. Let us imagine these subspaces as situated at the vertices of the 
square lattice on the torus of the sizes I x m (which will mean the periodic boundary conditions 
in both discrete space variables). Let the subspaces be arranged in checkerboard fashion, as in 
Fig. 1.1, where the empty circles correspond to subspaces of the space Tit, while the filled circles 



correspond to those of the space 7i 2 . 

Let then the operators A, . . . , D be such that the image of each of the mentioned d-dimensional 
subspaces with respect to, say, operator A lies in the d-dimensional subspace of Ti.2 at which points 



the arrow marked "A" that links these two subspaces (Fig. 1.1). Analogously, the restrictions on 
B, C, D are depicted in Fig. [hi] (see also formula (3.11) for non-degenerate A, . . . , D). Thus, to 
each link of the lattice a d x d matrix is attached that is a block of one of the "large" ma trices 
A,...,D. Let us shade half of the squares of the lattice in a checkerboard way, as in Fig. 1.1. One 
can verify that the evolution of the system may be described as follows. 

At the first step, each of the four d x d matrices that correspond to the arrows surrounding 
each shaded square is transformed into a matrix expressed through just these four matrices. This 
goes according to the following formulae, in which the d x d blocks are somewhat freely denoted 
by the same letters A, . . . , D as the "large" matrices: 

(3.1) 

(3.2) 
(3.3) 
(3.4) 



A - 




BD- X C) 


B - 




AC~ l D) 


C - 


-(C- 


DB- X A) 


D - 




-CA- X B) 



G 



However, the formulae (3.1-3.4) apply equally to the "large" matrices. 



After the transformation (3.1-3.4), all the arrows reverse, and at the second step the non-shaded 



squares are engaged in the same way according to the same formulae (3.1-3.4). Then everything 



is repeated. Thus, the evolution is of hyperbolic nature: each local perturbation spreads not faster 
than one unit of length per unit of time. 

Let us clarify the symmetries of vacuum curves and divisors Dx in this "reduced" model. Let 
us introduce two integer-valued coordinates £, r\ for the vertices of the lattice, so that £ increases by 
1 in passing from a vertex one step to the right, and r\ increases by 1 in passing one step upwards. 
£ and rj are defined modulo I and m respectively. A d-dimcntional subspace of Tii or Ti.2 will be 
denoted Ti.^ v if it corresponds to a vertex with coordinates £, rj. Consider a linear transformation 



in spaces Hi and Ti.2 consisting in multiplying the vectors of each subspace H^ v by i 
a fixed primitive root of the Z-th degree of unity: 



u>\ being 



1. 



This corresponds to the following transformation of the operators A, . . . , D (from now on we speak 
of each of these operators "as a whole", not of their blocks): 



A^uj x A, B^B, C^C, D-hj^D. 



(3.5) 



Consider also another linear transformation in Tii and Ti.2, consisting in multiplying the vectors of 
each subspace 7i^ v by tu^ , u>2 being a fixed primitive root of the m-th degree of unity: 

wj* = 1. 

This corresponds to the following transformation: 

A -> A, B — > cj 2 B, C -> w^C, D —> D. 
The vacuum curve of the operator L, which is given by equation ( |2.4| ) 

P(u, v) = det(uA + B- uvC - vD) = 0, 



(3.6) 



must be invariant under the transformations ( |3.5| ), (3.6). This leads to the invariance of the 
polynomial P(u,v) with respect to the following transformations g\ and <?2 : 



gi(u,v) = (uj 1 u,uj 1 1 v), 



(3.7) 
(3.8) 



This invariance, then, leads to th e fo llo win g statement: only those coefficients cijk are non-zero in 
the vacuum curve equation (see (2.2), (2.3)) for the "reduced" model, for which 



j — k = 0(mod I), 
j + k = 0(mod m) 



(3.9) 



As for the divisor Dx , let us recall that it consists of such points in the curve V in which vanishes 



the last coordinate of the vector X (see 42 ) , the latter being an eigenvector of the matrix ( p.6[ ) 
with zero eigenvalue: 

(uA + B -uvC -vD)X(u,v) = 0. (3.10) 

This immed iate l y le ads to the conclusion: the divisor Dx is invariant with respect to the trans- 
formations (fTiJ U). 

Under some additional condition, the inv erse st atem ent also holds: if the curve T and divisor 
Dx are invariant under the transformations (3.7), (3.8), then the corresponding i-operator comes 



from a "reduced" model described in this section. For instance, this is true if 1/2 and m/2 are 
relatively prime numbers. If these numbers are not relatively prime, some conditions are to be 
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imposed on the divisor -Co- 
here. 

Thus, let an operator L 



To avoid going into details of this latter case, let us not consider it 



.4 
C 



B 
D 



be given, A, . . . ,D being nxn matrices, n = (lm/2)d, I 



and m even, and 1/2 and to/2 being relatively prime. Let the vacuum curve T of the operator 
L and the divisor Dx be invariant under the action of the group Q generated by its elements 
3i)ff2 (3.7, 3.8), u>i and loi being primitive roots of degrees I and to of unity. Then the linear space 
in which operators A,...,D act decomposes into a direct sum of lm/2 d-dimensional subspaces 
7i.£ri, £ and r\ being integers modulo I and to respectively and such that £ + r\ is an even number, 
and the following equalities between the images of these subspaces hold (in a "generic" case of 
non-degenerate A, . . . , D); 

AH^-i,n+i = BH^n = C7~t£,ri+2 = DH^ + i. n +i- (3-11) 

The equalities ( |3.11 ) mean exactly that one is in the situation of Fig. 1.1. 

Let us prove the above statements. First, the natural projection from the curve T to its factor 
Y/Q has no branch points (here the fact that 1/2 and to/2 are relatively prime is used to demonstrate 
that ramification does not occur when aor« equals zero or infinity). Thus, the n-dimensional 
linear space of meromorphic functions x(z) — x(u, v) whose pole divisor is Dx decomposes into 
a direct sum of subspaces of equal dimensions corresponding to the characters of (commutative) 
group Q. Each of these subspaces consists of functions x{z) satisfying relations 



the character x^r) being a scalar factor 



where 



£a rjb 



9192- 



The equality g\ g™^ 2 = 1 means that £ + 77 must be an even number. 

The components of the vector X(z) are exactly the functions x{z). In an appropriate basis, 
d components correspond to each character ■ Let us denote H^ v the set of vectors with other 
components equal to zero. Now, the equalities (3.11) are to be proved to end this section. 

Consider the decomposition of vector X(u, v) into a sum 



X(u,v) 



where X^, v 6 l~t£, v - Then 



Consider the sum 



E 

geg 



Xtv(9^)9{{uA + B-uvC- vD)X(u, v)} = 



(3.12) 



(which is equal to zero because of ( 3.1C )). The action of g upon the braces in (3.12) means that 
each u and v in the braces is transformed according to (3.7), (3. 
and v changes into x—1,-1 
the order of group Q, 



i.e. u changes into xi,-i{g) u i 
(g)v. The equality (3.12) gives thus, after cancelling a factor equal to 



!iM^_i., |+ i(u,t)) + BX^ iV (u,v) — uvCX^, v+ 2(u,v) — vX, 



(u,v)D = 0. 



(3.13) 



Let us set it = in (3.13). Then v can take n different values vj satisfying relation P(0,vj) = 
0. To these values vj correspond d linearly independent vectors X^ v (0,Vj), and also d vectors 
Xt+i. v +i{0,Vj). Thus, the equalities 

BX iv (0,Vj) = VjDX£ + i, v+ i(0,Vj) 



that result from ( 3.13 ) give 

Analogously, one can as well obtain the rest of equalities ( [3.11 ) . 
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-vd i ua 






H 


b- 1 

) wa X 


r ^ 

—uvc 


7 —1 


—uvc 


— H 


—uvc 




) — — — H 

— UVC 

^vd_j 


H 




T 





Figure 1.2. The ways along these arrows are connected with both the vacuum curve and the dimer 
model 



4. Connection to dimer model 

As has been demonstrated, the integrals of motion of the dynamical system of Section [j] and its 
reductions (if the even degrees of the transformation (1.1) are considered) are the coefficients a,jk 



of the vacuum curve (2.3). These coefficients are determined up to a common factor, so they may 



be divided by aoo- As one can see, the resulting coefficients are those of the polynomial 

det(l + uAB- 1 - vDB- 1 - uvCB- 1 ). (4.1) 



In other words, the determinant (4.1) is an integral of motion for any u, v. 

Let us turn now to the model from section |], that is to the model in 2+1-dimensional discrete 
space-time with periodic boundary conditions, and let the dimension d of the linear space corre- 
sponding to each vertex be equal to 1. Each of the "small" matrices A,B,C,D corresponding to 
the links will then be a single (depending on the link) number a, 6, c or d. It is well known that 
the determinant of any N x N matrix is a sum of its matrix elements products corresponding in a 
certain way to the permutations of N objects, while each permutation decomposes into a product 
of the cyclic ones. In our situation, the cyclic permutations correspond to the non-sclfintcrsecting 



closed paths (contours) going along the arrows of the following diagram (Fig. 1.2) (thus, general 
permutations correspond to the sets of non- intersecting paths). To each closed path corresponds 
the product of the weights ua, —uvc, —vd, b~ x on its links, and, to get right signs for the terms 
of which the determinant ( |4.1[ ) is made up, one should add a minus sign to each such product 
containing an even number of the factors b^ 1 . 

Remark 4.1. Another way to obtain right signs is: to multiply each b by — 1 and then multiply 
each product corresponding to a closed path (and containing any number ofb's) by —1. 



It turns out that the determinant (4.1) is connected with the statistical sum of the well known 
dimer model f39f . Let us define the correspondence between the sets of paths and the dimer config- 
urations as follows. Let the empty set of paths correspond to the "standard" dimer configuration, 



the dimers being placed on the "S-links" (Fig. 1.2). For a non-empty set of paths, let us change 



the standard configuration along all the paths, replacing each dimer by a free link and vice versa. 
One can verify that this is a bijective correspondence. 

The statistical sum being considered, let the weights —b (not 6 _1 ) correspond to the "S-links", 
while to the other links correspond the unchanged weights ua, —vd, —uvc. Then one can see that 
the statistical sum, if multiplied by Jlover all links (~ (^ et us can the result the normalized 



statistical sum), consists of the same terms as the determinant (4.1), up to different signs of some 
of them. Let us emphasize that the dimer model is, of course, inhomogeneous: the weights a, b, c, d 
are different for different links. 

Let us study these signs in detail. Note that the conditions of non-intersecting and non- 
sclfintersecting impose strong restrictions on the possible path configurations. Every closed path 
on the torus is homologically equivalent to a linear combination with integer coefficients of two 
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„ L 



Figure 1.3. The standard dimer configuration 



basis cycles a and b whose intersection number is 1 (I use the boldface font for cycles, because the 
letters a,b . . . are already in use). If the torus is cut along a closed non-selfintersecting path c not 
equivalent to zero, the result will be homeomorphic to the lateral surface of a cylinder (this follows, 
e.g., from |68| , chapter 1, section 3). Then the contour d going along a generatrix of the cylinder 
in a properly chosen direction has the intersection number 1 with the contour c. The intersection 
number being bilinear and integer-valued, we find that if the contour c is homologically equivalent 
to a sum la. + mb, then I and m cannot have common divisors (not equal to ±1). Thus, the 
following lemma is valid. 

Lemma 4.1. Every closed non-selhntersecting path on the torus is homologically equivalent to a 
linear combination of the basis paths a and b with relatively prime integer coefficients. 

Now let us pass to the case of several contours on the torus. If they do not intersect, their 
intersection numbers equal (of course) and thus their homological classes must be proportional 
to one another. This together with Lemma 4.1 leads to the following lemma. 



Lemma 4.2. Several closed non-intersecting and non-selhntersecting paths going along the arrows 
on the torus, as in Fig. [Q| are necessarily all homologically equivalent to one another. 



If two paths are homologically equivalent, then the terms of the same degrees in u and v 
correspond to them (one can see in Fig. 1.2 that the different ways round an "elementary square" 
yield the same degrees of u and v). Let the basis paths a and b yield the terms proportional 
to x = u otl v^ 1 ,y — u a2 v^ 2 correspondingly (with the factors of proportionality not depending on 



u,v). According to Lemma 4.1, the determinant (4.1) and the statistical sum of the dimer model 
are polynomials in x, y. The following lemma sums up this section. 



Lemma 4.3. Let f(x,y) and s(x,y) be the determinant ( 4.1 ) and the normalized statistical sum 
of the dimer model considered as functions of x and y. Then 



s(x,y) = ^(-f{x,y) + f{-x,y) + f{x,-y) + f(-x,-y)), 

/(>> v) = \ (- s ( x ' y) + s (~ x > v) + s ( x > -v) + s (~ x > -y)) ■ 

Proof. If the normalized statistical sum consists of the terms 



(4.2) 
(4.3) 



then the determinant consists of the same terms multiplied by 



(-1) 



number of contours 



= r_i)g- c -d-(i> fe ) — (_iy k +i 



(here Remark 4.1 and Lemmas 4.1 and 4.2 are used). This means that the signs of all the terms 
must be changed except where both numbers j and k are even. This is exactly what the formu- 



lae (1.2, 4.3) do. The lemma is proved. 
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5. Equation of motion for "physical" variables 



For the dynamical system connected with the dimer model, to each link of the square lattice 
corresponds a complex number denoted a, b, c or d, as was explained in the beginning of Section |J 
Recall that the dependence of those numbers on a link of the lattice is implied, but usually not 



indicated explicitely, not to overload the formulae and figures such as Fig. 1.2. We study the 
evolution of the dynamical system up to the gauge transformations ( |l.3| ). It is easily seen that, 
for the dimer model, those transformations can be described as follows. Let us call the elementary 
gauge transformation corresponding to a given vertex of the lattice and a non-zero complex number 
A the following: take the numbers a, b, c and d corresponding to the four incoming and outgoing 
links for the given vertex, and multiply them all by A. The general gauge transformation for the 
dimer model will be a composition of any number of elementary transformations (which, of course, 
commute with each other). 



Remark 5.1. There exists also a discrete group of transformations of the form ( p.4 ) compatible 
with the reduction of the system to the dimer model which consists of translations of the lattice: 
G = H^ 1 in ( |l.3p IS m this case the operator of translation by an integer vector (£, rj), with £ + n 
even. Here we, however, do not consider such transformations. 

In this section, we will introduce, instead of a,b,c,d, new variables £1 which can be called 
"physical" in the sense that fl does not change under (general) gauge transformations. We will 
write out the "equation of motion" for fi, which turns out to be a 2+1-dimensional variant of the 
Toda lattice in discrete time. The remarkable and not a priori evident property of this equation is 
the completely equal status of the three coordinates: spatial £, r\ and the time r. We will consider 
also the questions of how to come back from f2's to the variables a, b, c, d, and whether it is possible 
to express the statistical sum in terms of Q. 

So, let a number f2 correspond to each two-dimensional cell (square) of the lattice, which is a 
"multiplicative circulation" composed from the numbers a, b, c, d on the links bordering the cell as 
follows: 

n = a ea b eb c e "d Ed , 

where 

£<i = —£b = —£c — £d = ±1, 

the sign "+" corresponding to non-shaded squares on Fig. |l.l| , and the sign "— " corresponding 
to shaded squares. In other words, let us move around the square anticlockwise, and take the 
number on each link in the degree one or minus one if we pass that link in or against its direction 
correspondingly. 

Remark 5.2. This construction is described naturally in terms of cohomology theory. Namely, 
the torus T 2 with the square lattice on it and the given orientation (direction ) for all links of the 
lattice is a CW-complex (163/, chapter 1, §4). The "Geld" a, 6, c, d is a 1-cochain with coefficients 
in the multiplicative group C* of non-zero numbers (we assume that in the "generic" situation, 
which we prefer to consider, no one of the numbers a, b, c, d equals zero), while its "rotor" O is the 
coboundary of that cochain. This cohomological interpretation will be of use for us later, when we 
will consider the reconstruction of the "held" a, b, c, d from Q. 

Let us examine the change of circulations fl in one step of evolution. Here one must consider 
separately the squares of the form c 



a 
d 



b (i.e. shaded ones in Fig. 1.1) and the squares b 



d 
a 



(non-shaded in Fig. 1.1). The former we will call the active squares, the latter — the non-active 
ones. 

The circulation around an active square is 

fl = a^cd^b. (5.1) 



11 



The transformation (3.1 



3.4), which in our case can be written as the transformation 



a — > d {ad — be) , b^—c(ad — bc) , , , 

c — > — b (ad — bc)^ 1 , d^a(ad-bc)^ 1 

with the simultaneous reversing of all the arrows, changes the circulation in an obvious way as 

ft^ft- 1 . (5.3) 



Convert now the formulae (5.2) for active squares to the following form which we need for 

(5.4) 



examining what happens in the neighboring non-active squares: 

a -» a- x (l- ft)" 1 , 6^ &-!(!- ft- 1 )- 1 



c — > c 



Hi -ft- 1 )- 1 , d^d-^l-Q)- 1 . 



This done, consider the formulae (5.4) not for one active square, but for four active squares adjacent 
to the chosen non-active one. To be exact, let us consider the formula a — > a- 1 (l — ft) -1 for the 
active square having the common link a with the chosen non-active square, and similarly treat 
the transformation formulae for b, c and d, so that a,b,c,d are now the numbers corresponding 
to the four links surrounding the non-active square, while the four ft's in ( |5.4| ) become different, 
namely the circulations around the four active squares. Now we see that the circulation around 
the non-active square itself 

ft = ft(£,r?) = ac^db- 1 , 

£, rj being the integer-valued coordinates of, say, south-western vertex of this square, are trans- 
formed by one step of evolution as 

ft(£, 77) - ft(£, r,) ■ (1 - ft(£, 77 - 1)) • (1 - ft(£, 77 + 1)) x 

x (1 - ft- 1 ^ - l^))- 1 ■ (1 - ft- x (£ + M))" 1 . (5.5) 

After the step of evolution, the active squares become non-active and vice versa. Thus, with 
the proper choice of the origin for the integer-valued time coordinate r, the active squares have 
the odd sum £ + 77 + r of their three coordinates. As for the non-active squares, they can be totally 
eliminated from consideration when deriving the "equation of motion" for ft because, according to 



(5.3), the circulation around a non-active square is simply the inverse of the circulation around the 
same square in the previous moment r, when it was active. Writing out explicitely the dependence 
of ft on the time, we get from ( [5.3| ) and ( |5.5| ) the following "equation of motion" containing only 
the circulations around active squares: 

ft(e,r ? ,r-l)ft(e,r;,r+l) = 

_ (l-n(g,q-l,T))(l-n(g,77+l,T)) 
(1 - ft-HC - 1, 77, r)) (1 - ft-!(£ + 1, 77, rj) ' 

£ + rj + t even. 



Note the following compact form of equation (5.6). Introduce the "discrete pscudo-Laplacians' 



A Ti £ and A,,^ acting on functions F on the cubic lattice, by the formulae 

(A T , S F) (£, 77, r) = F(£, 77, r - 1) + F(£, 77, r + 1)- 

-F(£-1,7/,t)-F(£ + 1,77,t); (5.7) 

(A„, e F) (£, 77, r) = F(£, 77 - 1, r) + F(£, 77 + 1, r)- 

-F(£-1,77,t)-F(£ + 1,77,t). (5.8) 



By means of these operators, ( p.6|) is rewritten simply as 

A^lnft = A,, iC ln(ft - 1). (5.9) 
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It is seen from (5.9) that the spatial coordinate n and the time r in our model are present at an 
equal status. Namely, the equation (5.9) obviously does not change under simultaneous changes 



T <-> 77, 



ci <-► 1 -n. 



(5.10) 



The equal rights between the two spatial axes must exist, of course, too. It is described, as one 
can easily verify, by the substitutions 



V, 



ci <-> rr 1 . 



(5.11) 



It is implied, of course, that the interchanges of coordinate axes in ( |5.10 > 5.1l| ) act only on the 
subscripts of pseudo-Laplacians of the type ( |5.7| , 5.8). Those interchanges generate the symmetric 
group 53 acting on the set (£,77, r), and the corresponding transformations of CI generate a group 
of linear- fractional transformations isomorphic to S3. 

Remark 5.3. The equation or (fjj can be obtained by a simple transformation from the 



"discrete Toda Geld" equation from p2\j (see als o \5% p(\ pl\l), at least if we do not take into account 
the boundary conditions. That discrete held in \62j depends on an integer-valued coordinate k and 
two more coordinates u and v, each taking values with steps h, and is written as f k (u,v). The 
field equation has the form (formula (8) from jj&sj/) 

exp(/ fe (w + h,v + h) - f k (u + h,v)- fk(u,v + h) + f k (u,v)) = 
= (l + h 2 cxp(/ fe+ i(u + h,v) - f k {u,v + h))^ x 

x(l + h 2 exp(f k ('u + h,v)-f k - 1 {u,v + h)) y j. (5.12) 



Divide ( 5.12 ) by the equation obtained from ( 5.12 ) by the change 
k — ► k — 1, u — > u — h, v — > v + h, 
and introduce a new variable 

9k{u,v) = -h 2 e-Kp( fk(u + h,v) - f k -i (u, v + h)) . 
We then get the equation 

g k (u,v + h)g k (u-h,v) _ (l - flfe+i(u, v)) (l - g k -i(u - h, v + h)) 



9k(u, v) g k (u-h,v + h) (l - g k (u, v)) (l-g k (u-h,v + h)) 



Finally, setting 



£ = -h 1 u + h 1 v, 



h- 1 



v, 



to{€,v>T) = gk(u,v), 

we get for VI equation ( |5.6| ). 

Thus, the dynamical system connected with the inhomogeneous dimer model is a variant of the 
Toda lattice. 

Let us now pay attention to the question to what degree can the "field" a, b, c, d be reconstructed 



from its "rotor" Q. First of all, this can be done only if Q is indeed a coboundary (see Remark 5.2), 
i.e. if the corresponding to CI cohomology class in the group H 2 (T 2 , C*) vanishes. The mentioned 
cohomology group, according to the book [ p8| , is isomorphic to C*. Constructively speaking, all 
this looks as follows: two "fields" CI belong to the same cohomology class if and only if the products 
of their values over all two-dimensional cells are equal (recall that the group C* is multiplicative), 
and the field CI is a coboundary if this product equals 1. 

Assuming that this last condition is satisfied, we can find a cochain a, b, c, d whose coboundary 
is Cl. It is natural to want to determine that cochain to within the (general) gauge transformations 
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introduced in the beginning of this section. It is easy to see that, in cohomological terms, the 
gauge transformations are changes of the 1-cochain a, 6, c, d by coboundaries of 0-cochains. From 
its coboundary, however, fl can be restored only to within a cocycle. To find a, 6, c, d up to a 
coboundary, one must fix two more parameters because, according to p8| , the group i? 1 (T 2 , C*) is 
isomorphic to C* © C*. Those two parameters are nothing else than the values of cochain a, 6, c, d 
on two (arbitrarily chosen) basis cycles on the torus, i.e. the products of numbers corresponding 
to the links, taken in degrees ±1, over all links entering a basis cycle. Here, as well as earlier, the 
sign + is taken if the direction of our way around the cycle coincides with the direction of a link, 
and the sign — is taken if it does not. 



Comparing this situation to the representation of the motion integrals (4.1) as a sum over the 
sets of closed paths on the torus given in Section ||, we see that the determinant (4T) and, conse- 
quently, the statistical sum are not determined uniquely by a field Q. They are determined only 
up to a two-parameter transformation which, as one can see from the considerations of Section |], 
can be interpreted as a renormalization of the variables u, v: 

u — > const! • it, v — > const2 ■ v. 



6. The discrete analog of Lax pair and a generalization of 
the dynamical system 

Now let us return from the reduction of Section [| to general matrices L = ^ ^ ^J^J ' us 
consider the evolution described in Section |l| from another viewpoint. Denote 



This means that 



A B 
C D 



( A 








\ B 


d ; 


(c 


S) 



(6.1) 



It follows from the equality (6.1) that 

AA + CC = BB + DD, (6.2) 
AB + CD = 0, 
BA + DC = 0. 

These three equations are equivalent to the fact that the following equality holds for any complex 
u: 

-(A- uB)~ 1 {C - uD) = (uA + B)(uC + (6.3) 
Vice versa, from ( |6.3| ) follows 

LL ~ { F 
F being equal to both sides of (|6.2|), i.e. 

It is clear that with any choice of F the matrix L belongs to the same equivalence class. The 
formula (3.3) defines the same evolution in the space of these classes as it was in Section ^, with 



the agreement that the operators without a tilde correspond to the moment of time r, while those 
with a tilde corre spo nd to the moment t + 1. 

The formula ( |6.3| ) suggests the following generalization. Let, from now on, A(u) and B(u) be 
matrices depending polynomially on u: 

A(u)^A + A 1 u + --- + A mA u mA , (6.4) 
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B(u) = B + B x u + ■ ■ ■ + B„ lB u mB . (6.5) 

We will look for matrices A(u), B(u) — the matrix polynomials of the same degrees via and in 
u — that satisfy, for any u, the equation 

= A^Biu)- 1 . (6.6) 



The relation (6J3) provides what is called a discrete analog of the Lax L,v4-pair. Let us forget 
for a moment that we are already using the letters L and A for other purposes, and remind that 
the Lax L, A pair is a pair of operators depending on the time and other parameters and having 
usually some special form, the evolution of operator L in the case of discrete time being described 
by the formula 

L(t + 1) = A{t)L(t)A(t)- 1 . (6.7) 



On the other hand, let us rewrite (returning to our notations) ( |6 -6[ ) in the following form: 

Biu)- 1 !^) = A^Biuj^A^Aiu)- 1 . (6.8) 

Comparing (6.7) and (6.8), we see that, in our case, B(u)~ 1 A(u) plays the role of Lax L-operator, 
while A(u) — the role of Lax A-operator. 

Let v be an eigenvalue of both sides of ( |6.6[ ). Let Y(u,v) be the corresponding eigenvector 
normalized, as in Section ^, so that its last coordinate equals unity, and let X{u,v) be the vector 
proportional to B(u)~ 1 Y(u, v) and normalized in the same way. One can verify that this may be 
described by the following formula {h(u, v) being a scalar factor): 

( ^] ) X(u, v) = h(u, ») ( J ) ® Y(u, v), (6.9) 



which is in obvious analogy to ( |2.5| ). The divisor equivalence is 

mD u + D x ~D v + DY, (6-10) 

D u and D v being pole divisors of the functions u and v, m — max(rriA, tub)- 
For a given u, the eigenvalues v come from the equation 

P(u,v) = det(A(u) - vB(u)) = 0. 

It defines an algebraic curve T — "generalized vacuum curve". Let us calculate the genus g of the 
curve r. First, we need to know the number of branch points of the projection 



(u,v) — >u (6.11) 



of the curve T onto the complex plane. 
Consider P(u,v) as a polynomial in v. 



P{u,v) = a (u) +ai(u) H h a n (u)v n . (6.12) 

One can verify that a,j (u) has a degree 

deg aj(u) = (n - j)m A + j B - (6.13) 

From this one can deduce that the discriminant of P[u,v) considered as a polynomial in v is a 
polynomial of degree 

b = (rriA + ms) n (n — 1) 



in u. The mapping ( 3.11 ) being n-sheeted and the number of branch points equalling b, one obtains 
from the Riemann — Hurwitz formula that 

9 = (n-l)(^±^n-l). (6.14) 
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So, the following construction has been described. Given two polynomial matrix functions A(u) 
and B(u), one considers the meromorphic matrix function A(u)i?(w)~ 1 (or else B(u)~ 1 A(u)), and 
from this function the algebro-geometrical objects arise: the generalized vacuum curve T and 
the linear equivalence class of the pole divisor Dy (or, respectively, Dx) of the eigenvectors of 
the mentioned meromorphic matrix function. Instead of the pair (A(u), B(u)), it is sufficient to 
indicate its equivalence class with respect to gauge transformations 

A(u) -> GA{u)H, B(u) -> GB{u)H; (6.15) 

instead of the function A(u)B(u)~ 1 , its equivalence class with respect to transformations 

A{u)B{u)- 1 -> GA{u)B{u)~ 1 G~ 1 

will suffice. Then it turns out that the correspondence between such equivalence classes (either of 
the pairs (A(u), B(u)) or the functions A(u)i?(w) _1 ) and the abovementioned algebro-geometrical 
objects is a birational isomorphism, the divisors D x and Dy being of degree g + n — 1, as in 
Section ||. 

The easiest way to show this is to start from a given curve T defined by the equation 

n (n-j)m.A+jm B 

P(u,v)=^2 a jk v j u k = Q 

j=0 k=0 



(compare with ( 6.12 , 6.13 )) and a divisor Dx in it of degree g + n—1. The number of coefficients 
a,jk minus one common factor equals 

(n + 1) f!M+H n + 1 )_i, (6 .i6) 



The linear equivalence class of divisor Dx is defined, as is known, by g parameters. Adding up 
the expressions ( |6.16| ) and (6.14), one gets the total of 



(mA + m B )n 2 + 1 (6.17) 



parameters. 



Then, the gauge equivalence class of the pair (A(u),B(u)) is constructed out of re latio n (6.9). 



To give more details, one must at first choose a divisor Dy satisfying the equivalence ( 6.10 ). Then 
the poles and zeros of the function h(u, v) are determined. For X(u, v) and Y(u, v) one must take 
columns consisting each of n linearly independent meromorphic functions with corresponding pole 
divisors. The arbitrariness in these constructions leads exactly to the fact that A(u) and B(u) are 
determined up to a transformation ( |6.15 ). 



The pair (A(u), B(u)), up to a scalar common factor, is determined by (yriA + raj + 2)n — 1 



parameters (see (6.4, p.5| )). In taking the gauge equivalence class, the number of parameters is 
reduced by 2(n 2 — 1). The result is again (6.17). This means that, indeed, to a generic pair 



(A(u), B{u)) corresponds a divisor Dx of degree g + n—1 and the correspondence 

(gauge equivalence class of the pair {A{u),B{u))) < — > (r, class of D x ) 

is a birational isomorphism. 

Now let us recall that Y(u, v) was defined as an eigenvector of the operator ^4(u)i?(it)~ 1 , while 



X(u,v), as is easily seen, is an eigenvector of B(u) 1 A(u). The relation (3.6) means that for the 



pair (A{u),B{u)) its vector X(u,v) is nothing else than Y(u,v), i.e. the equivalence holds 

D x ~D x + {mD u -D v ). (6.18) 

Now, assuming that if a quantity without a tilde corresponds to the moment of time r then 
that with a tilde corresponds to r + 1 , one comes to a conclusion that to the adding of unity to 



the time corresponds a constant shift ( 3.18] ) in the Jacobian of the curve T. Thus, the dynamics 
of the system in this section, as well as in Section ^|, linearizes. 
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Chapter 2 

Dynamical system connected with 
the transposition of three matrices 



7. Definition of dynamical system 



Let 



A 




(7.1) 



be a block matrix acting in the linear space of m + n + r-dimensional complex column vectors, 
so that, for example, A, F and K are square matrices of sizes mxm, n x n, r x r respectively. 
Consider the problem of factorization of matrix A into a product of the form 



A = A X A 2 A 3 , 



(7.2) 



where 



At 



At Bx 
Ci D 1 
1 



A 2 



A 2 B 2 
10 

C 2 D 2 



10 

A 3 = | A 3 B 3 
C 3 D 3 



(7.3) 



The factorization ([7^) may be seen as a generalization of the factorization of an orthogonal rotation 
in the 3-dimensional space into rotations through the "Euler angles" . However, at the moment we 
consider no orthogonality conditions. 



One can obtain from the factorization (7.2) other factorizations of the same kind by using the 
following transformation of the triple Ai,A 2 ,A 3 : 






f 


1 









M 2 












M 3 


1 





N 







1 





) 








M, 1 , 





A 3 



(7.4) 



M\, M 2 and M 3 being arbitrary non-degenerate matrices of proper sizes. 



Lemma 7.1. For a generic matrix A factorization (7.2) is (if exists) unique to within the trans- 
formations (7.4). 
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Proof. First, let us specify that matrix A will be called generic with respect to this lemma if it 



a) is non-degenerate and b) allows a factorization (7.2) with matrices A 2 and T> 2 non-degenerate. 
Let 

A = AAA (7-5) 



be another factorization of the same kind. From (7.2) and (7.5) one finds 



A = AAA : 



3 > 



(7.6) 



where 

A( = (A 1 )- 1 A 1 , A^A^AV 1 - (7.7) 

Let us denote the blocks in the dashed matrices by the same letters A%,. .., D3 as in equalities ( |7.3[ ) 
with proper number of dashes added to them. The relation (|7.6|) is rewritten as 



(7. 




KM 


A'{B 2 C^ + B'iA 


A!lB 2 D'i + B'{B'i 


C'{A 2 


C'{B 2 C% + D'lA 


C'{B 2 Dl + D'lB'i 


c 2 


D 2 C'i 


D 2 D'i 



B' 2 


D' 2 

From here, one obtains at once the equalities 

C' 2 = C 2 , C[=0, C"J = 0. 
Taking this into account, one finds from the block in 2nd row and 2nd column that 

D»A» = 1. 

Thus, D'l and A are non-degenerate. Now the blocks just above the main diagonal yield 



B'l = 0, 



So, the matrices A'{ an d A ( [7-7|) are block-diagonal. In is easy to see that this means exactly 
that A[, A' 2 , A are obtained from Ai, A 2l A3 by the transformation (7^). The lemma is proved. 

Now let us construct, starting from the block matrix A, new matrix B by following means: 
factorize A into the product (|7.2|) and set 



B = A 3 A 2 A 1 . 



(7.9) 



From the above considerations it is seen that the matrix B is determined to within the transfor- 
mations 








AI. 






Mo 



(7.10) 



Let us call such transformations, as applied to the block matrices here, the gauge transformations. 
The following simple but important observation is valid: if matrix A itself undergoes a gauge 



transformation, this in no way affects the set of matrices B obtained from formula (7.L). 

It will be shown in Section |l0| that factorization ( |7.2| ) does exist for a generic matrix A. This 
factorization will be constructed by means of algebraic geometry. Taking this into account, we 
are ready now t o d efine the dynamic system that we are go ing t o examine. Let M. be the set of 
block matrices (7T) taken to within gauge transformations ( 7.1C ), or, using stric ter language, the 
set of equivalence classes of such matrices with respect to transformations ( 7.10 ). The set M. will 
be our "phase space" . Then, the birational mapping / is defined on the set M. that brings into 
correspondence with a matrix A, factorized into the product (7.3), the matrix B factorized into 
the product ( |7.9| ). Let us now bring into consideration the "discrete time" r taking integer values 
and say that to the trasition from time t to time t + k corresponds the mapping 
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8. Connection with Clifford algebras and "twisted" Yang- 
Baxter equation 

The dynamical system described in Section |t] deals with linear operators acting in the direct sum of 
three vector spaces. Now we will show that one can pass on from the direct sum to a tensor product 
of (other) spaces, and connect our dynamical system with the "twisted" , or generalized, Yang- 
Baxter equation (about which see [|50[ §2|). The arising solutions of that equation belong to 



the "free fermion" type. This will be also seen in Section 13, where we, basically, will also arrive at 
a "twisted" Yang-Baxter equation. In the present section, we consider the most general situation, 
and use algebraic means, while in Section [l] we will consider some reduction of our model, and 
will use some topological considerations as well. 

Let us begin with the constructing, out of a given square matrix R of size (n\ +712) x ("4 +^2), 
of the operator 1Z acting in the tensor product of spaces £\ and £2 having dimensions 2 ni and 2™ 2 
correspondingly. 

Consider the Clifford algebra generated by fermionic creation and annihilation operators a^, 
where 1 < j < N, i.e. the algebra with relations 

a t a t + a t a t = a 7 a k + a k a J = Oi 

_ + + _ r ( 8 - x ) 

a j a k + a k a j = 6 3k, 

Sjk being the Kronecker delta symbol. Such algebra can be realized by operators in a 2 N - 
dimensional space, and than there exists in that space a non-zero vector ("vacuum") 51 such 
that 

ajfl — for all j. 

Let us have constructed two such algebras, with N = n\ and n 2 , acting in spaces £\ and £ 2 - 
We are going to explain how to obtain out of them, in the spirit of |33| , one "large" algebra acting 
in £ 1 (gi £2 ■ 

Consider the particle number operator 

"1 

A/i = 5>+aJ (8.2) 

i=i 

acting in £\. The eigenvalues of this operator are integers, so one can introduce the operator 

Zx = {-l)*K (8.3) 
Then let us define the operators 6^ in £ 1 ® £ 2 by formulae 

bf =af®l if 1 <j <m, (8.4) 

bf =Zt® af_ ni if m + 1 < j < rix + n 2 . (8.5) 

(it is understood that the left factor in the tensor products acts in £\, while the right one — in £2). 
A simple check shows that the operators (8.4- §7^) satisfy the same commutation relations ( |8.l[ ), 
and the vacuum is now the tensor product of vacua in £\ and £ 2 ■ 

Let us now associate to a matrix R the automorphism ip of the Clifford algebra acting in £ 1 <g> £ 2 
which acts on the creation and annihilation operators in the following way: 

711+712 

k=l 



n 1 +n 2 



fc=i 
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According to the theory of Clifford algebras (as it can be found, e.g., in pH]), <p is an inner 
automorphism, i.e. it can be represented in the form 



(p(A) = 7tA7T 



(8.6) 



A being an arbitrary element of the algebra, 1Z s Aut £\ ® £2 ■ To conclude the construction of 
operator 1Z, it remains to define, similarly to (B.2), the particle number operator 



AA 



ri2 



in the space £2, and to set 



ft = •£•(-!) 



(8.7) 



Theorem 8.1. Let the matrices A\, A2, A3 be as in Section [?j (formulae (7.c )), and let there be 
given also the "primed" matrices A'i, A' 2 , A' 3 of similar form, e.g. 

' A[ B[ 0' 
A[ = \ C{ D[ 
1 



etc. Let the equality 



A1A2A3 = A'zA^A'i. 



hold. We will associate to the matrix Ai or, to be exact, to its "nontrivial" part R = 



A t 
C'i 



Di 



the operator 1Z = 1Z\2 acting in £\ (g> £2 according to the construction described by formulae (8.1 
8.7), setting n\ — m, 712 = n. In the same way, let us associate to the matrix A[ the operator 

n' = n' 12 . 

Consider also a space £3 of dimension 2 r and, replacing £ 2 by £3 and n% by r in the construction 



ft8.l\-p. 7y, associate to the matrices A 2 and A' 2 the operators C = £13 and £ = C 13 , correspondingly 
(the letter C plays the role oflZ from (8j_ ), the subscripts of C and £ denote the numbers of spaces 
£j in whose tensor product an operator acts). Finally, to the matrices A3 and A' 3 we will associate 
similarly the operators M = M13 and M' = M.\ 
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Under these conditions, the "twisted" Yahg-Baxter equation holds: 

KCM = M'C'K', 



(8.9) 



where we imply, as usual, that each operator is multiplied tensorly by the unity operator in the 
lacking space, e.g. C is multiplied by 1 € Aut £ 2 ■ 



The equality ( p.8[ ) must be interpreted as the decomposition (7.2) of the operator A obtained 
at the previous step of evolution from Section [?] as the product 

A = A'sAA 



of type ( f7.9D - Thus 
a reformulation of ( 



can be regarded as a description of a step of that evolution, while (8.9) is 



in which the direct sums of vector spaces are replaced by tensor products. 



Proof of Theorem 8.1. To begin with, introduce operators TZ, C, Ai, Ai', £ , 1Z' as follows. 



Generalize the construction (8.4-3.5) for the tensor product £ 1 <8> £ 2 (8) £3 of three spaces, joining 



the creation and annihilation operators acting in them in a single algebra by multiplying them, if 

}]) hold. We get 



necessary, tensorly by operators of type (S.3), so that the upper relations 



.1 



cf = Z X 



of ® 1 ® 1 if 1 < j < m, 
(g) o^- ® 1 if m + 1 < j < m + n, 
' Z% <S> a j if m + n + l<j<m + n + r. 



(8.10) 
(8.11) 
(8.12) 



20 



The operators TZ, . . . , TZ' in space £\ <£> £2 <£> £3 are determined by the following conditions: 1) each 
of them doesn't change the vacuum Sli ® ® ^3 and 2) consider a column of height m + n + r 
made of the annihilation operators cj (8.1C- 8~L2] ). The action of each of the operators TZ, . . . , TZ' 
on that column by conjugation of every element, e.g. 

is equivalent to the action of the corresponding matrix Ai , A2 , A3 , A s , A 2 or A[ on the whole 
column, e.g. 

ci \ / Ci 

-> Ai 



Such operators TZ, 
Besides, the equality 



, TZ' , according to the theory of Clifford algebras ]56j], exist and are unique. 



TZCM = M'C'TZ', 



.13) 



holds, because both sides of 
operator column, namely: 

ci 



13|) preserve the vacuum and act identically on the annihilation 



A1A2A3 



Cl 



m+n+r - 



A 3 A 2 A 1 



Cl 



The following lemma shows how the operators in ( 3.15 ) (they don't have, ge nera lly, the form 
required in Yang-Baxter equation) are connected with the operators in relation ( 3.S) that we are 
proving. More exactly, we will connect the operators "with hats" with the operators "with tildes" 
which, we remind, were defined each in its own tensor product of two spaces. 



Lemma 8.1. 



TZ' = TZ, 



K = TZ, M — M, M' = M', 
while C and £' are connected with C and £ as follows: 

where (/) denotes the prime or its absence, and N1N2 is, strictly speaking, A/i <& A/2 ® 1. 
Proof of Lemma 



(8.14) 
(8.15) 



is based on the fact that LHS's and RHS's of flS.14| ) and ( p.!5| ) define 
(by conjugations) the same automorphisms of the Cliffor d algebra acting in £\ ® £2® £3- This is 
check e d by direc t calculations for each of the equalities ( S.14 , 8.15 ) and for each set of o perat ors 



10, B.ll, 8.12) separately. We will remark only that the role of factors (— 1) 



A/i A/i 



seen from equalities of type 

(_l)A/-xA/- 2 ( Zi ® ± ® 1) = (1 ® a^= ® 1) (-l)^ 2 , 



A§ is 



(8.16) 



which demonstrate how an operator, in this case, Z\, disappears or appears after a conjugation 
with (-l)^ 2 . 



To conclude the proof of Theorem 8.1, let us write out the equality (3.13) in terms of operators 
TZ,...,TZ' defined by equalities of type (8.7), using (8.14) and (8.15): 



72,£(_^)A/i(A/i+A/'3)^ / (^_ ^\A/iA/"3 _ 

= M'(— c 1 {— i) A/ "i( jV ' 2 + A ''3)7^'(— 1)^1-^2 



.17) 



All the degrees of minus one in LHS and RHS of ( p. 17 ) can be conveyed through other factors 
to the right using the fact that M. commutes with M% + A/3 (i.e. "preserves the total number of 
particles of second and third kinds" ) and with A/i , similarly £ commutes with N\ + A3 and with 
A/2, and TZ' commutes with A/i + Af 2 and with A/3. This done, the factor (_l)MAr 2 +MA4+AA 2 A4 ) 
arising in both LHS and RHS, is canceled, and we come to (I8.9D . The theorem is proved. 
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9. Invariant algebraic curve of matrix A and some divisors 
in it 

The dynamical system of Section turns out to be completely integrable. To be exact, an invariant 
curve r can be constructed out of matrix A, together with a divisor T> on it. In terms of these 
algebro-geometrical objects, the evolution is as follows: T does not change, while V — more precisely, 
its linear equivalence class — depends linearly on the discrete time r. 

Let us start from the definition of the curve T. The word "invariant" in this definition will be 
justified in Section [l0|. 



Definition 9.1. The invariant curve T of the operator A of the form (7A) is an algebraic curve in 
CP 1 x CP 1 x CP 1 (i.e. in the space of three complex variables u, v, w, each allowed also to take 
value oo) given by equations 

det(A- I vl n I) = 0, (9.1) 



ul m 











Vl n 











wl r 


V 


= uw. 





(9.2) 

Here a subscript of each 1 means the size of corresponding unity martix, while denotes 
rectangular zero matrices of different sizes. Strictly speaking, equations (9.1, 9.2) define the "finite 
part" of the curve T, the whole curve T being its closure in Zariski topology. 

_ fx 

The equality ( ) means that a column vector X = Y ] exists, with X, Y, Z column vectors 

V z 

of dimensions m, n and r, correspondingly, such that 







( uX 




- 


vY 






\ wZ 



A\ Y = vY . (9.3) 
\Z J \wZ J 

Such vectors X form a one-dimensional holomorfic bundle over T. 

The next lemma shows the structure of the zero and pole divisors of functions u, v, w. For these 
divisors, the notations (u)o, (u)oo etc. are used. 

Lemma 9.1. There exist such effective divisors (i.e. hnite sets of points) T>i, . . . , T>q in the curve 
T that 

{u)oo=V 1 +V 2 , {v) 00 =V x +V z , {w) 00 =Vz + X> 4 , 

(tt) = V A + £> 6 , («)„ = £> 5 + V 6 , Ho =V 5 +V 2 . { ' ' 

V3 and V§ are of degree m, T> 2 and V4 are of degree n, T>\ and Vq are of degree r. Generally, 
all points included in divisors T>\ , . . . , 2?g are different from each other. 




Proof. Co nsider, e.g., the case u = 0, w 7^ 0, w ^ 00. Then, according to (9.2), v — 0. The 
equality (9.3) turns into the following system: 

BY + CZ = 0, 
FY + GZ = 0, 
JY + KZ = wZ. 

One can express X and Y through Z (e.g., Y = —(F - DA- 1 B)- 1 (G - DA- 1 C)Z) and then 
substitute these expressions into the third one. One will come to an equation of the form 

KZ = wZ (9.5) 

which has r characretistic roots Wx, . . . , w r , different from each other in general case. This is how 



r points (0,0, u>i), . . . , (0,0, uv) of divisor Vq are obtained. The other divisors in (9.4) arise in a 
similar way. The lemma is proved. 



22 



The vector 




in (9.3) is determined up to a scalar factor which may depend on the point 



in the curve T. So, this vector can be normalized by setting its hrst coordinate identically equal to 
unity (cf. J42[). Y becomes then a meromorphic vector on T with a certain pole divisor T>. 

v z J 

However, X, Y and Z taken separetely satisfy stronger restrictions, as the following lemma shows. 
In the lemma, (/) denotes the divisor of a function / (zeros enter with the + sign, poles with the 
— sign, as usual). 



Lemma 9.2. The column vector X consists of functions f such that 

(/) +V- («)«, > 0; 

the column vector Y consists of functions f such that 

(J)+V-(v) oo >0; 

the column vector Z consisits of functions f such that 

(f)+V-(w) oc >0. 



(9.6) 
(9.7) 
(9.8) 

Proof. One can see immediately from the formula ( |9.3| ) that the vector uX entering into R.H.S. 

X 

cannot grow faster than the vector | Y | in L.H.S. in such points where u = oo. This is exactly 

Z 



what the inequality ( |9.6| ) states. The inequalities (9.7) and ( |9.8| ) are proved similarly. 

Now the time has come to make it sure that the curve T, for a generic matrix A, is a smooth 
irreducible curve. One may wish also to calculate its genus in some simple way. To do that, we 
are now going to examine a relatively simple particular case of the matrix A, although at the time 
"generic" enough to make sure that such its features as genus and the degree of divisors are the 
same for matrices in some Zariski neighborhood. 

Thus, let all the matrix elements of A equal zero except the ones lying, hrst, in the main 
diagonal and, second, in the "broken" diagonal parallel to the main one (for these latter matrix 
elements, the difference between the numbers of a column and a row must be some constant modulo 
m + n + r). The elements in the main diagonal will be denoted as a%, .. . , a m , fx, . . .,/„, fci, .. . ,k r ] 
and let the elements in the broken diagonal be all equal to the same complex number s: 





/ ai 




S 




\ 


A = 


s 


h 

fn 


S 










h 










S 




k r J 



(9.9) 



It does not matter through which blocks exactly the "broken" diagonal passes. 

For the hnite u, w, the curve T now examined is given by equation (resulting from the substi- 
tution of Q and (|9^) into ([□])) 



F(uw) = H (a a - u) J] (f fi -««,)• J] (*7 " w ) ± gm+n+f 

a=l 0=1 7=1 



0. 



(9.10) 
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As is known, in singular points 

dF/du = 0, 



dF/dw = 0. ( 9 ' n ) 

The system ( 3.11 ) has a finite number of solutions, and, changing s in (|9.10| ), one can make these 
solutions not to lie in the curve T, which thus will be free of singularities for finite u, w. It is an 
easy exercise to show that there are no singularities when u or w is infinite as well. 

Retu rning now to general matrices A and curves T, let us note that it cannot be that the 
system ( |9.11 ) or its substitute in the neighborhood of infinite u or w possesses solutions in the 



curve r in general case and does not possess them in a particular case. Thus, the smoothness of 
r for a generic A is clear. As for irreducibility, to prove it let us examine the natural projection 
of r onto the Rimann sphere CP 1 of the variable u. This projection is an (n + r)-sheet cover, 
and if T consisted of two or more components, the sheets of the cover would split into groups 
belonging to each component. To prove that it is not so in the general case, it is enough to 
present an example where it is not so. To do this, take A of the form (B^9) and, moreover, put 



/i = • • • = f n = k± = ■ ■ ■ = k r = 0. Equation (9.10) then becomes 



Y[ {a a -u)± s m + n + r = 0. 



a=l 



Let a\ ^ and not coinside with other a a . Then in a neighborhood of the point (u, w) — (aj., oo) 
the variable w behaves, up to a nonzero factor, like 



_1 - (ai - u) 



l/(n+r) 



From here one sees that all the mentioned n + r sheets belong to a single component, i.e. the 
irreducibility of T is proved. 

Now let us denote the number of branch points of the cover T ^ C 3 u as b. Then the genus 
of the curve T, according to the Riemann — Hurwitz formula, is 

g = l-n-r + ^. (9.12) 

Our next aim is to express b and g through m, n and r. 

X 

Lemma 9.3. The degree of the vector [ Y | pole divisor T> is m + b/ 2. 

Z 



Proof. Write out the equation (9.3) "explicitly' 



AX + BY + CZ = uX, 

DX + FY + GZ = vY, } (9.13) 
HX + JY + KZ = wZ. 

Expressing X through Y and Z by means of the first of these equations and substituting into the 
rest, one finds: 



it 



1 (D(u- A)- 1 B + F) u- 1 (D(u~A)- 1 C + G)\(Y\_ (Y 



H(u-A)- l B + J H{u-A)- 1 C + K )\Z ~ \Z 



w _ . (9.14) 



To a generic u correspond n + r different w = w%, . . . , w n + r , and the corresponding n + r vectors 
Y 



Z 



are linearly independent as eigenvectors of the matrix in L.H.S. of (9.14). An easy check 



shows that in the "suspicious", from the standpoint of equation (9.14), points u = 0, oo, and also 
in points where det(w — A) = 0, there exist as well n + r linearly independent vectors ( „ ] 



-solutions of limit cases of the system (9.14). 
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Consider a determinant 

Y(wi) ... Y(w n+r ) 



(9.15) 



Z{W\) . . . Z(w n+r ) 

Given u, it changes its sign under odd permutations of w's. 

However, d 2 is a function of u only. From the above one sees that d 2 (u) vanish in branch points 
where to a given u correspond less than n + r values of v or w. This yields b zeros of the function 



d 2 (u). There is, however, one more cause for this function to vanish. According to Lemmas 9.2 
and |9.l| , Y and Z must vanish as a whole in the points of divisor T>^ — the common part of divisors 
{v)oo and (u>)oo- The degree of divisor P3 is m, so we get 2m more zeros of d 2 {u). 

The number of function d 2 {u) poles equals the number of its zeros, i.e. b + 2m. Thus, the 
Y 
Z 



mcromorphic vector ( _ I has 6/2 + m poles, as desired. It remains just to note that vector X 



Y 



in accord with formula (9.6), has no poles that ( ~ ) doesn't have, and on the other hand doesn't 



vanish in points of divisor P3. The lemma is proved. 

Let us turn again to matrices A of the form ( |9.9[ ) . For such matrix, it is easy to find the vector 

Y = X in a given point (u, v, w) £ T. Let us assume that the "broken" diagonal is the one 

Z ) 

adjacent to the main diagonal, so that there is only one letter s in the lower left corner. Then the 
following holds for the vector X coordinates: 



(01 


- u)X x 4 


-sX 2 -- 


= 0, 




- u)X 2 4 


-SX 3 : 


= 0, 




- u)X m 


+ sYi 


= 


(A 


- v)Y l + 


sY 2 = 


0, 




- v)Y n 4 


sZy = 


= 0, 


(h 


- w)Z 1 4 


- sZ 2 = 


= 0, 




- w)Z n - 


hsX! 


= 0. 



From here the ratios between vector X coordinates are readily seen. Assuming the normalization 
condition X\ = 1, one finds out that X has the poles a) of the order m in n points (u,v,w) = 
(oo,//3,0) and b) of the order m + n in r points (u,v,w) — (00, 00, fcy). In all, X possesses 



thus mn + mr + nr poles, taking their multiplicities into account. Recalling Lemma 9.3 and 
formula ( 9.12| ), one can now find the genus g of the curve as well. As a matrix A of the form fl9.9Q 



is "generic enough" , the results on the degree of divisor T> of the vector X and genus g of the curve 
apply also to curves corresponding to generic matrices A. Let us formulate them as the following 
lemma. 

Lemma 9.4. For a generic matrix A, the genus of the curve T is 

g = mn + mr + nr — m — n — r+l, (9.16) 



while the degree of divisor V of the meromorphic vector X = Y is 




mn + mr + nr — g + m + n + r — 1. (9-17) 

Thus, in this section we have constructed, for a given matrix A, an algebraic curve T and a 
bundle of vectors X over it, and calculated the genus g of the curve and the degree of the bundle 



(i.e. the divisor V degree). As a helpful tool, a matrix A of special simple form (9.E) was used 
which, from many viewpoints, was "generic enough" . In Section |l^ we will study how these objects 
behave under evolution introduced in Section 
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10. Evolution in terms of divisors 



In this section it is shown, at first, that there exists a one-to-on e cor respondence (more precisely, 
a birational isomorphism) between the set of block matrices A (7T) taken up to gauge transfor- 
mations (7.4), and the set of pairs (an algebraic curve, a linear equivalence class of divisors on it) 
of a certain kind. This correspondence has, in essence, been constructed in Section ^, and here 
are some missing details. Then, it is explained which divisors and why correspond to the factors 
A\,A2 and A3 in (7/2) taken separately. Finally, it is demonstrated that to the matrix B (7.9) 
obtained from A by reversing the order of its factors, the same curve F corresponds, but the divisor 
undergoes some constant shift. Thus, the motion linearizes in the Jacobian of curve T. Let us 
proceed to a detailed consideration. 

Equations (9.1, 9.2) define, for a block matrix A, an algebraic curve Y. Those equations can 
obviously be written as 



m n r 



a ijkU i V j W k = 0, 

j=0 j=0 fe=0 



(10.1) 



X 

Besides, a linear bundle over Y has been constructed in Section |— the bundle of vectors Y 

Z 



(9.3). That means that the divisor T> of the bundle is determined, up to linear equivalence, whose 
degree is g + m + n + r — 1, g being the curve's genus ( 9.17 ). Gauge transformations (7A) do not 
change a pair (r, class of divisor T> ). 



Now let us show how to construct the matrix A starting from coefficients aijk of the curve (10.1) 
(arbitrary complex numbers in "general position") and a divisor T> of degree g + m + n + r — 1. 
Note that genus g of the curve Y defined by formulae (10.1) without any (a priori) connection 
with block matrices is given by the same formula (9.16). This can be seen, e.g., by starting again 
from the "simple" curve ( 3. 1C ) of Section ^| whose genus is known. Define now the meromorphic 
column vectors X, Y and Z, guided by Lemma |9.2| : for components of vector X, take m linearly 
independent meromorphic functions on Y satisfying relation ( |9.6| ), and for Y and Z take, similarly, 
n functions satisfying ( |9.7D and r functions satisfying ( |9.8| ). 

Note also that Lemma 9.1 about divisors (u)oo, (v)ooi (w)oo entering in formulae (|9~6-3. j) re- 



mains valid for curves defined by an "abstract" system (10.1), which is immediately seen on sub- 
stituting zero or infinity for u, v, or w in (10.1). 

It is clear now that relation ( |9.3j ) determines unambiguously the matrix A (cf. a similar con- 
struction in paper [jL?| ) . Another choice of linearly indep endent functions for components of X, Y 
and Z leads, of course, to a gauge transformation ( 7.10[ ). The vectors X, Y, Z change under it to 
M\X, M2Y, M3Z. On the other hand, if divisor V is changed to another divisor belonging to the 
same linear equivalence class, the vectors X, Y and Z are just multiplied by a scalar meromorphic 
function h having zeros in the points of the first divisor and poles — in the points of the second one. 

Hence, the following theorem has been proved. 



Theorem 10.1. If block matrices A and A of the form (7.1) are connected by a gauge transfor- 
mation 



A 









(10.2) 



M 3 



then they have a common invariant curve Y given by equations of th e fo rm ( 10.1 ), and the holo- 
morphic bundles of the vectors X corresponding to the m b y formula ( |9.3| ) are isomorphic. 

Conversely, if two matrices A and A of the form ( 7A_ ) have the sam e invariant curve Y and 
the corresponding vector X and X bundles are isomorphic, then ( 10.2 ) holds. Being properly 
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Figure 2.4. Factorization of matrix A and the divisors 



normalized (say, by a condition that the first coordinate identically equals unity), X and X become 
meromorphic vectors with linear equivalent pole divisors T> and T> of degree g + m + n + r — 1, and 



X 




h(u, v, w) 



(10.3) 



where h is a scalar meromorphic function whose divisor (h) = (h)o — (h)oo satisfies equality 

(h) =V-V. 



Examine now each multiplier in factorization (7.2) separately. Lemma 7.1 shows that factor- 
ization (7.2), if exists, is unique to within the transformations (7.4). Let us demonstrate how to 
construct this factorization by algebro-geometrical means. 

Consider the following figure (Fig. 2.4). 

The meaning of the numbers standing near the edges in this figure is as follows: if those numbers 
are jk, then the meromorphic vector corresponding to the edge consists of such functions / whose 
zero and pole divisor (/) satisfies inequality 

(/)+£>- v 3 -V k >0. 



Those inequalities must be in agreement with Lemmas 9.1 and 9.2. In particular, the matrix A3 
will be defined by equality (notations of formulae (|7.3j) are used), 



A 3 B3 

C3 As 



Y 
Z 



Y' 
Z' 



(10.4) 



where the meromorphic vector Y consists of functions / such that 

(/) +V - V 1 - V z > 

(formulae fl9.7p and (9.4)); Z of functions such that 

(/) + V - V 3 - V 4 > 
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(formulae (9.8) and (9.4)); Y' and Z' consist by definition of such linearly independent functions 
that 

(/) + V - V 3 - V 6 > 

for Y', and 

(/) + V - V 2 - V 3 > 



for Z' . It is easy to see that ( 10.4 ) is a correct definition for matrix A 3 , because the components 

Y \ , / Y' , 

form a basis in the space of meromorphic functions / 



of each of the vectors 
such that 

Next, let 



and 



Z 1 



(f)+V-V 3 > 0. 



A 2 B 2 
C 2 V 2 

A x Bi 

where X' consists of functions / such that 

(/) + 



X 
Z' 

X' 
Y' 



X' 

wZ 

uX 
vY 



(10.5) 
(10.6) 



V - V 2 - P fi > 0. 



It is shown in much the same way as above that equalities ( |10.5[ ) and ( 10. 6| ) do correctly define 
the matrices A 2 and A\. What remains is to check the validity of equality ( |7.2| ) for A\, A 2 and 
^3 given by these definitions. To do this, observe that ( 1 . 4| — |ToT6| ) together yield 



A1A2A3 







( "X \ 


I) 




vY 






I, wZ / 



(10.7) 



The equality (7.2) follows from comparing (10.7) with (|9.3| 



Note that the arbitrariness in choosing X',Y' and Z' corresponds, of course, to transforma- 



tions (7.4) 



Now l et u s pass to matrix £>, a product of the same three factors in th e inv erse order. The 
and ( 10.4 - 10^ ) together yield (if one multiplies both sides of ( 10. 5| ) by u, and both 



formulae 



sides of (10. 4|) by v 



B 







1 uX 




- 


vY 


uZ 1 ) 




\ vZ 1 



(10.8) 



Compare the divisors of meromorphic vectors in L.H.S.'s of (10. S) and ( |9.3| ). An easy calculation 
shows that 

T> x > -V x = V Y , -V Y = V (uZI) -V z =V x - V 6 . 

One sees hence that the same curve T corresponds to the operator B as to the operator A, while 
the divisor V changes to V + T>\ — T> 6 . 

Thus, in this section the name "invariant" has been justified for the curve T: it has been shown 
not to change under the evolution of Section [?]. At the same time, it was demonstrated how to 
construct the factorization (7.2). Finally, it was shown that the evolution is described in algebro- 
geometrical terms as a linear, with respect to discrete time, change of (the linear equivalence class 
of) divisor T>: it changes by 2?i — Vq per each unit of time. 
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11. Algebro- geometrical objects in the orthogonal and sym- 
plectic cases 

In this section, the algebro-geometrical devices developed earlier are supplemented with necessary 
means for studying the evolution of orthogonal and symplectic matrices A. It turns out (Section |l2| ) , 
that the condition of orthogonality or symplecticity is an admissible, i.e. compatible with the 
evolution, reduction of the dynamical system defined in Section ^. We will start, however, with 
some consideration for generic matrices A, namely, looking at what happens with the invariant 
curve r and the bundle of vectors X (see formula ([)!])) under two operations: inversion A — » A^ 1 
and transposing A — > A T . 

As for the inversion, here virtually everything is seen from formula ( |9.3| ), especially if one 
rewrites it as 

A' 1 I vY ] = I Y | . (11.1) 
\wZ ) \Z ) 

It is clear that the curve T, basically, does not change; the roles of column vectors X, Y, Z are 
played by uX,vY and wZ, while the roles of meromorphic functions u,v,w on T are played by 



u 



■ v 



and w . The exact formulation is given in the following lemma. 



Lemma 11.1. If to a matrix A corresponds 
v = uw (see (fO 



9.2 ) or (10.1)), and a divisor V in it (see Lemma 9 
it), then to the matrix A~ l corresponds the curve 1^ 



curve r given by equations P(u,v,w) — and 
and the paragraph before 



given by equations P(u 



" 1 ) = 



and 



v 



uw, 



and the divisor T>( ^ in it that is the image of T> under the natural isomorphism 

r^r ( - 1} : (u,v,w) i-» (m -1 ,!! -1 ,™ -1 ). (11.2) 
Proof. It remains to remind that T> is the pole divisor common for the meromorphic column 



vectors in LHS and RHS of (9.3) or (11.1). Hence the validity of the statement in lemma about 
the divisor Z?( -1 ) is clear, and thus the lemma is proved. Let us add, however, that, for the curve 
r( -1 ), the roles of divisors £>i, Vq, T>2, 2?4, "Da and T>^ (let them be listed in this order) are played 
by the images of divisors Vq, T>\, £>4, T>2, T>^ and D3. This is due to the fact that the role of 
divisor (u)^ is played by the image of (u)o, the role of (u)q is played by (u)^, and sim ilarly for 
the functions v and w. The connection between those divisors is given by formulae (3.4). 

Consider now the transposing A — > A T . The determinant is not changed under transposing, 
thus it follows from (9.1) that the invariant curve of matrix A T is the same T. Somewhat more 
complicated is to find column vectors X, Y, Z, satisfying the equation 



A 1 







( uX 




- 


vY 






\ wZ 



(11.3) 



We will obtain them here using the direct constructive method from the author's work |pl[ . Note 
however that we might act also in another way, adapting to our case the idea of divisor duality and 
of the scalar product of meromorphic functions on an algebraic curve associated to this duality, 
ascending to the work |34[ and described also in more recent works 65, ^6|, |67fl . We will pay 



and described also in more recent works 65, 
tribute to this elegant idea in Section |l8| in the framework of the "local" approach to the evolution 
of orthogonal matrices, while returning now to the methods of ]8T 
Transpose (11.3): 



X T Y T Z T 



A = (uX T vY T 



i.e. by the column vector 



Multiply both sides of ( |11.4| ) from the right (scalariy) by 

3D, and use the formula (9.3) to exclude A from the LHS. One gets, after collecting 



appearing m 
similar terms, 




(11.4) 



(u - u)X T X + {v- v)Y T Y +{w- w)Z T Z = 0. 



(11.5) 
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In formula (11.5), (u, v, w) and (it, v, w) are so far arbitrary points in the curve T (one and the 
same, as was explained above). Set now u = u. To a given (generic) it correspond n + r different 



u>'s (recall the formulae (10.1)). Choose now w and w different: w ^ w. Then (11.5), on being 

w 

multiplied by -, yields 



vY T Y + wZ x Z = 0, where v 



uw. 



Let w be fixed, and w take other n + r — 1 possible values which we will denote W2, ■ ■ ■ , w n + r . 

vY , 

are proportional to the cofactors of the determinant first 



Then the components of vector 



wZ 



column entries 

A(it, w) 

Choose the normalization as follows 



Y{u,UW,w) Y(u,UW2,W2) 
Z(u,UW,w) Z(u,UW2,W2) 



Y(u,uw n+ri w n+r ) 

Z(u,UW n +r,W n +r) 



vY 
wZ 



A(it, w 



(vector of the 
first column entries 
cofactors of A(it,-u;) 



(11.6) 



(11.7) 



We have already met the determinant A(it, w) under the name of d (formula Here, 
however, we consider it not as a two- valued function of a point (it, uw, w) G L rather than a function 
of it. This leads to a large multiplication of this determinant zeroes, some of them being also zeroes 
of all the minors — components of the vector from (11.7), while others being not. The zeroes of 



A(it, w) are of interest for us, of course, as candidates for being poles of the vector ( ^ ^ ) . Note 



that, although the enumerator and denominator in the RHS of (11.7) are two-valued functions, 
the resulting fraction is obviously single- valued. 

Consider first the "branch points", i.e. the points where either a) Wj = w for some j, or 
b) Wj = Wk for some unequal j and k. There are b points of the first type (see formula ( |9.12| )), 
and those points form the ramification divisor Z? ra m, while the points of the second type are of 
no interest for us, because all the minors in such points have, too, zeroes of the same character 



(namely as the square root of a local parameter in a curve), so that zeroes in the RHS of (11.7) 
cancel one another. Next, we know that in some points the columns of determinant A(it, w) must 
vanish as a whole. Such points, too, are divided into a) those where the first column vanishes and 
b) others. The first are again m, and they form the divisor T>3, while the second are again of no 
interest, because the minors also vanish in them. 

( vY \ 

Other candidates for being the poles of the vector I ~ g J might, be the poles of the enumerator 



in the RHS of ( |l 1.7 ). However, in a point where a column of minors has a pole, the corresponding 
column of A(it, w) has a pole, too. Recalling also the linear independence of columns of A(u,w) 
(see the proof of Lemma 9.3), we get a pole of A (it, w) that cancels the pole of minors. 

' vY ' 



Thus, the poles of the vector I J can onr y be the points of the divisor 

£>ram + 2?3 • 

We must consider also points where all coordinates of this vector vanish. Of all candidates the 
only survivers here are the first order poles of A(it, w). There are ran + mr + nr such poles, and 
they form the divisor denoted in this paper as V. Other candidates in zeroes of all coordinates 
vY , 

are sifted away by arguments like those used above when searching for the 



of the vector 



poles. 



wZ 



Summarize these considerations in the following lemma. 
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Lemma 11.2. The vector 
relation 



vY 
wZ 



from (11.7) consists of functions f whose divisors obey the 



(/) + ©re 



V 3 -V>0. 



It remains to clarify the following. Those functions, in principle, may obey some stronger 
constraints, such as 

(/) ~t~ ©unknown ^ 

0, Where £> un known < ©ram + V 3 -V. 

It turns out that they do not. This is seen from the degrees of divisors: 

deg (Z?ram + ©3 — ©) = b + m — (mti + mr + nr) = g + n + r — 1 
(recall ( |9.12 ) and ( |9.16 )), and this is exactly what we need. 



X 

Denote now as © the divisor of singularities common for the vectors | Y 

Z 



and 



(we 



mean here that p oles ente r with the plus sign!). As, for the divisor T> and vectors X, Y, Z, the 
same relations as (9.6-?.8) hold, one can derive that T> is obtained from the divisor of singularities 

/ gy\ i — i 

of the vector I _ ^ , by adding T> 5 (for illustration, Fig. 2.4 again can be used). Thus, we got the 
\wZ J 

following: under the transposing of a matrix A, the divisor T> changes to 

© = © ram + ©3+© 5 -©- (11.8) 

Here 2? ram is the ramification divisor with regard to the variable u. In the following theorem, this 
result is formulated in more symmetric and elegant form. 



Theorem 11.1. Let to a matrix A correspond an algebraic curve V and a divisor V in it (as was 
described in Sections ^ and [h]), while to the transposed matrix A T — (the same curve and) a 
divisor T>. Then 

V + V~ V can +V 1 +V 2 +V 3 +V 4 + V 5 + X> 6 , (11.9) 



where T> can is a canonical divisor in T, and divisors T>\, . . . , T>q are defined in Lemma 9.1 . 

Proof. As is known [ )78| , canonical divisor is the divisor of any differential in T. Take the 
differential du and find its divisor (du) = (du)o — (du) oo- H is not difficult to understand that the 
zeroes of the form du are situated exactly in the branch points of T considered as a covering over 
the Riemann sphere of variable u: 

(du) = ©ram, (11.10) 



while the poles of du coincide with the poles of function u, but are of order 2: 

(du) oo = 2{u) 0O = 2{V 1 +V 2 ), 



(11.11) 



according to Lemma 9.1. It follows from (11.10) and (11.11) that, for any canonical divisor V c 

T> Tam ~V can + 2(Vi+V 2 ). 



For m ore symmetry, recall that T>\ + T>2 ~ (u) z 
( 11.9 ). The theorem is proved. 



(u)o ~ 2?4 + 2?6, after which ( 11.8 ) turnes into 



Theorem 11.2. If a block matrix A of the form ( |7J[ ) is such that A T = MA~ 1 M~ 1 , with M 
having the form 

Mi 
M 2 
M 3 



M 



(11.12) 
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(in other words, A T and A 1 are gauge equivalent, cf. the formula ( 7.1(fy ), then the invariant curve 
r of matrix A possesses the involution 



I : (u,v,w) < — ► (u ,v ,w 1 ), 
while the divisor T> and its image D 1 under involution I obey the equivalence 

V + V 1 ~£> can + Pi + ... +V 6 . 



(11.13) 



(11.14) 



Converse ly, if a matrix A is such that its curve T possesses involution ( 11.18 ), and the divisor 
equivalence ( 11.14 ) holds, then A T and A -1 are gauge equivalent. 



Proof follows at once from Lemma 11.1 and Theorems 11.1 and 10.1 



Further properties of the block-diagonal matrix M. (11.12) providing the gauge equivalence 
between A T and A^ 1 come out during some algebro-geometrical examination, which we will now 
perform. These properties determine, roughly speaking, whether matrix A is "in essence" sym- 
plectic or orthogonal. Calculate M, using meromorphic vectors 



X 




and 



X 




of matrices A T and A 1 respectively, connected, according to Theorem 10.1 , by the relation 

X(u, v, w) — A4X(u, v, w)h(u, v, w) (11.15) 



(cf. formulae (10.2, 10.3)). To be exact, we will now deal with blocks M 2 and M3, so that the 
already obtain formulae (11.6) and (11.7) will be of use. 

For a given generic complex number u, there exist n + r points in the curve T with coordinates 
(u,uwi,wi), . . . ,(u,w n + r ,w n + r ). To avoid bulky formulas, we will omit the middle coordinate, 
and write those points as (u, wi), . . . , (u, w n + r ). One sees from ( 11.6| , 11.7) that Y and Z are given 
in all those points at once by the formula (where, of course, v\ — uw\ and so on) 



/ v\Y_{u,w\) 
\wiZ(u,wi) 



v n+r Y_(u,w n+r ) 
w n+r Z(u,w n+r ) 



Y(u,wx) 
Z(u,wi) 



Y(u,w n+r ) 
Z(u,w n+r ) 



For the vectors Y and Z corresponding to the inverse matrix, the relations 
vY(u, w) = y(u _1 , if -1 ), wZ(u, w) = Z(u~ x , w^ 1 ). 



(11.16) 



(11.17) 



hold. This follows from (11.1), if we replace there (u,v,w) by (u ,v , w ) (the replacement 
applies, of course, to the arguments of vector functions X, Y, Z implied in (11.1) as well). 
From the formula (11.15) follows (as before, we omit v in the triples (u,v,w)) 



( vY(u,w) 
\ wZ(u, w) 



M 2 OU vY(u,w) \ h( s. 
M 3 ){ W Z(u,w))- k{u > w) - 



(11.18) 



Giving w all n + r possible values, we obtain n + r linearly independent columns in the LHS and 
RHS, which allows us to express M 2 and M 3 from flll.l8| ). It is convenient to write the result, 
with equalities (11.16) and ( |11.17| ) taken into account, in the following form: 



M 2 ~ x 



Mo 



Y(u-\w^) 
Z(u-\w^ 1 ) 



Y(u-\w-l r ) 



Z[u 



-1 ,.1 ) 



) w n+r 
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/ h(u,wi) 
\ 



h(u, w 



Y(u,wx) 
Z(u,wi) 



Y(u,W n+r 
Z(u,W n+r ) 



(11.19) 



n+r ) 



Recall that h performs the divisor equivalence between V, given by the formula ( 11. 8| ), and V 1 
-the image of V under involution /, in the sense that 

(h) = V 1 - V. 



Substituting here ( |ll.8|) , we find 



(h) = V + V 1 - V lam -V 3 - V 5 . 



(11.20) 



The divisor in the RHS of (11.20) is obviously invariant with respect to involution /, thus the 
divisor (h) of function h possesses the same property. If the curve T is irreducible, and this is 
exactly the fact in the general position, then h is determined by its divisor up to a constant factor. 
This means that under the involution I the function h is multiplied by a constant which, evidently, 
must equal ±1: 



h(u,w) = ±h(u ,w 1 ). 



(11.21) 



Return to the equality (11.19). It has a constant matrix, not depending on a point of the curve, 
in its LHS. Take the matrix transpose of that equality, and change (u,w) <-» (u _1 , w^ 1 ). We get: 



M 2 






M 3 



± 



M 2 






M 3 



where the sign coincides with that in ( 11.21 ). It is clear that, for the similar reasons, the equality 
Mi = ±Mj also holds. As a result, the following lemma is proved. 

Lemma 11.3. If a block matrix A has the property 

A T = MA^M- 1 , (11.22) 

where M. is a block diagonal matrix, then M. is symmetric or antisymmetric if the function h in 
the curve T with the zero and pole divisor (11. 2C) is even or odd with respect to the involution I, 
correspondingly. 

It follows from (|ll.22j) that to a gauge transformation 

A -> NAN^ 1 , Af = diag (Ni,N 2 , N 3 ), 

corresponds the transformation 



M -> (Af- 1 ) MAf- 



(11.23) 



of matrix A4. A symmetrical AA can be reduced by such a transformation to an identity matrix, 
and for an antisymmetrical one each diagonal block Mi can be reduced to the standard form 

I consisting of half-size blocks. The relation ( [11.22; ) means, after this transformation, 

the orthogonality of A in the first case and its symplecticity — in the second case. 
Let us summarize these considerations in the following theorem. 



Theorem 11.3. In the notations and under the assumptions of Theorem 11. Z , a matrix A with an 
irreducible curve T is gauge equivalent to an orthogonal or symplectic matrix, if the function h with 
the zero and pole divisor ( 11.2(\ j is even or odd, correspondingly, with respect to the involution I 

Thus, in this section we studied the algebro-geometrical objects in a specific case of orthogonal 
or symplectic matrices A. Now it would be not very hard work to show that those objects retain 
their specific form under the evolution, thus proving that a matrix A retains its property to be 
gauge equivalent to an orthogonal or symplectic matrix. We will prefer, however, to go another 
way, considering the decomposition (7.2) in orthogonal and symplectic cases. 
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12. Matrix factorization in case of orthogonality or sym- 
plecticity, and conservation of those properties under 
evolution 

It is natural to expect that an orthogonal or symplectic matrix A can be factorized in a product 
of matrices of the same kind (orthogonal or symplectic) A\, A2, A3 (formulae (7.2, 7.3)). Indeed, 
the following lemma holds. 

Lemma 12.1. Under conditions of Lemma 11. c, for the matrix A such a factorization (7^ 7.3) 
exists that for each matrix Ai separately, i = 1,2,3, a relation like (11.22) is valid, i.e. 



Aj = MA7 1 M~ 1 . 



(12.1) 



Proof. Consider, to begin, an arbitrary factorization A = A1A2A3 of matrix A, not requiring 
that (12.1) hold. Then there are obvious factorizations 



and 



A T = AjAjAj 

*4 — 



(12.2) 



(12.3) 

One more factorization of matrix A T , besides (12.2), into a product of three matrices with zero 
and unity blocks in the same places, can be obtained from |T22| ) and ( )i~2~3| ) in the following way: 

A T = MA~ 1 M~ 1 = MA^A^A^M- 1 = 

= (MA^M- 1 ) (MA^M- 1 ) (MA^M- 1 ) . (12.4) 

Applying Lemma 7.1 to the two factorizations ( |12.2| ) and ( |12.4[ ), we deduce that the following 
relations, for some nondegenerate matrices F\, F2, F3 of proper sizes, must hold: 



10 

Aj = (MA3 X M- 1 ) I F 2 



F-i 



Ai 









(MA^M- 



Fo 




Ff 1 0' 
F 2 _1 I (MA^M- 1 ) 



(12.5) 



(12.6) 



(12.7) 











Our aim now is to find such a transformation of type (7.4) for matrices Ai, A2, A3 that the 
relations (12.1) hold. Consider at first a matrix 





A 3 = \ K 2 I A3 
K3 



(12. 



with some nondegenerate K2 and K3. Relations (12.8) and (12.5) together yield the following 
connection between Ai and A^ 1 : 



Ai = MAI 



This will turn into the required relation 















: 


K 2 







KJ 









Vo 













): 


F 2 









F3) 



(12.9) 



Al = MAn 1 M~ 1 



(12.10) 
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if we manage to obey, by a proper choice of K 2 and K 3 , the following equalities: 



K 2 M 2 1 K 2 T F 2 



Mo 



K 3 Mz l KjFz = Mz l . 



(12.11) 
(12.12) 



Rewrite ( 12.11 ) in the form 



K 2 M 2 x Kj = M 2 X F 2 X . (12.13) 
In the following Lemma 12.2] we will demonstrate that the RHS of ( 12.13j ) is symmetric or anti- 



symmetric in case M 2 is symmetric or antisymmetric correspondingly. Assuming, for a while, this 
fact without proof, we find that the equation ( 12.13 ) is always solvable with respect to K 2 — this 
is a simple consequence from the properties of quadratic forms and antisymmetric bilinear forms 
and their matrices, see, e.g., §§90 and 91 of the manual j7^] (one must take into account here also 
the nondegeneracy of matrices M^T 1 and F^ 1 ). 



Similarly, the equation (12.12) is solvable with respect to K 3 . Thus, we managed to obey the 
relation (12.1C). Then, solving the equation for K\, similar to (12.11) and (12.12), we see that the 
relations 



and 

also hold, for 
and 



= MA 2 T M~ 1 



Aj = MAjM~ 




Ai = A! 







•) 


A 2 \0 




V° 




1 
















1 








K7 



We don't need any more the initial matrices Ai wi thout tildes in the rest of this Proof. Rename 
the new matrices Ai into Ai, and with this Lemma 12.1 is proved. Recall, however, that its proof 
was based on the following lemma (and its analogs arising from changing the subscript 2 to 1 or 
3). 



Lemma 12.2. If, in the notations of Lemma 12.1 



M 2 = ±Mj 



then 

with the same sign as in ( 12.1^ ). 



{M^F, 1 ) = ± (M^Fz 1 ) , 



(12.14) 



(12.15) 



Proof. Take the equality (12.5) and apply to its both sides, first, the transposing, and second, 
the matrix inversion. We will get, in the LHS, A^ 1 , while in the RHS — some expression containing 
Aj , which we have no need to write down here. This done, express again Aj through A^ 1 by 
multiplying the LHS and RHS of the obtained equality by suitable matrices and the interchanging 
LHS and RHS. The result will be the following: 



AI 



MA3 1 





F T 
1 2 








F T 



M' 



Comparing (12.16) with (12.5), we find, in particular, that 

M 2 7 1 F 2 = F?MZ\ 



(12.16) 



(12.17) 
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It is easy to see that the equality ( [12.17 ) is equivalent to ( 12. 15j ) , and with this the proof of 
Lemma 12.2 is over, and the proof of Lemma 12.1 is thus complete, too. 



The following theorem, which summarizes the results of this section, immediately follows from 



Lemma 12.1 



Theorem 12.1. The property of a matrix A of the form (7A) to he gauge equi valent to an 
orthogonal or symplectic matrix and, consequently, to be factorable in a product (7J2, \7. <{ ) of 
matrices of the same kind, is conserved under the evolution described in Section 0. 

Proof follows at once from the definition of this evolution given in the end of Section 



3G 



Chapter 3 

Inhomogeneous 6- vertex model 



13. Dynamical systen connected with the 6-vertex model 
on the kagome lattice 

Consider now a reduction of the system denned in Section ^ which leads to a dynamical system 
in 2 + 1-dimcnsional fully discrete space-time. Let the linear space in which matrix A (7.1) acts 
have a basis enumerated by edges of a triangular lattice on the torus (Fig. |3.5[ ), so that those 

components of the vector X = Y corresponding to horisontal edges form the vector A, while 

those corresponding to oblique edges form the vector Y, and the ones corresponding to vertical 
edges form the vector Z. As the number of each of those three types of edges is the same, all the 
blocks in the matrix A have the same sizes: 



( = the number of lattice edges). 



(13.1) 



Impose the following "locality" condition on matrix A: let the vector corresponding to any given 
edge of the lattice be transformed under the action of A into a linear combination of just three 
vectors, corresponding to the edges coming upwards, to the right and northeastwards from the 
vertex that is the upper, right, or northeastern end of the considered "incoming" edge (Fig. 3.6). 
Thus, only those elements of matrix A are not zeros that correspond to "local" transitions of 
Fig. |] . 

The factorization of a "local" matrix A into the product (|7.2|) corresponds to each vertex 



represented by a small circle in Fig. 3.5 being converted into a triangle of the type shown in 



Fig. 2. 4, so that the lattice transforms into a kagome lattice (Fig. 3.7). The triangles arising from 
the vertices-circles are shaded in Fig. 3.7. 



One can easily see that the "locality" property of matrix A is preserved by a step of evolution, 
if, of course, the proper gauge is taken. The detailed description of the step of evolution from the 
"local" viewpoint and the description of gauge transformations preserving the "locality" of A are 
given in the beginning of Section Here we will concentrated on the values conserved under the 
evolution. 

Let us return to the triangle lattice of Figure 3.5. Express the "integral of motion" 



I(u,w) = det(l - A 





1 w~ 








(13.2) 



in terms of paths going along the edges of this lattice (I(u, w) is indeed an integral of motion 
with any it, w, because the equality I(u, w) — determines the invariant curve, and a possible 
multiplicative constant is fixed by the fact that the constant term in ( 13. 2| ) equals unity). 

As it known, the determinant of a martix is an alternating sum of its elements' products, 
each summand corresponding to some permutation of the matrix columns, while each permutation 
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Figure 3.5. The triangular lattice 




Figure 3.6. The vectors corresponding to "incoming" edges are transformed by A into linear 
combinations of those corresponding to "outgoing" edges 




Figure 3.7. The kagome lattice 
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factorizes into a product of cyclic ones. As applied to our matrix A, it means that the determi- 
nant (13.2) is a sum each term of which corresponds to a set of closed trajectories going along the 
arrows according to Fig. |3.6| (recall that the lattice is situated on the torus!). The trajectories of 
each given set can have intersections and self-intersections, but none of the edges may be passed 
through twice or more by one or several trajectories. 

To be exact, to each trajectory corresponds a product of entries of the matrix 



A 



-1 








u~ lr w 
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corresponding to transitions through a vertex to a neighboring edge according to Fig. 3.6, multiplied 
(the product as a whole) by (—1). To each set of trajectories (including, of course, the empty set) 
corresponds the product of the mentioned values corresponding to its trajectories. A direct check 
shows that all the minus signs, including that in formula ( |13.2| ), have been taken into account 
correctly. 

It is easy also to describe the determinant I(u, w) in terms of the kagome lattice obtained on 
factorizing the matrix A into the product ( |7.2| ). This description almost repeats two preceding 
paragraphs. Let us formulate it as the following lemma. 

Lemma 13.1. I(u,w) is a sum over sets of trajectories on the kagome lattice; the direction of 
motion is upwards, to the right, or northeastwards; none of the edges is passed through twice 

by trajectories of a given set; to the vertices of types , I , X , if a trajectory 

passes through them, correspond the factors equalling matrix elements of matrices A\, A2, A3 
respectively; besides, to each move to the right through a lattice period corresponds a factor it -1 , 
and to each move upwards — a factor w^ 1 (and both of them to a diagonal move); finally, to each 
set corresponds one more factor, (_l)(number of trajectories) _ 

Now let us link I(u, w) with the statistical sum of inhomogeneous 6- vertex model on the kagome 
lattice. Let each edge of the kagome lattice be able to take one of two states, which will be depicted 
below as either presence or absence of an arrow on the edge (the arrow will always be directed 
upwards, to the right, or northeastwards). A "Boltzmann weight" will correspond to each vertex 
as follows: if there are no arrows on the edges meeting at the vertex, the weight will be 1; if there is 
exactly one arrow coming into the vertex and exactly one going out of it, the weight will be equal 



to the corresponding matrix element of Ax, A2 or .A3 (e.g., to the vertex 



corresponds 



the weight equal to the matrix element of A2 that is responsible for a transition between the 
vector corresponding to the left edge, and the vector corresponding to the upper edge); if there 
are 2 incoming and 2 outgoing arrows, the weight is the difference between the products of weights 
corresponding to the intersecting and non-intersecting paths through the vertex: 



Weight 



= Weight 



Weight 



-Weight 



Weight 



(13.3) 



in the rest of cases the weight is zero. 

A weight will also correspond to each edge of the kagome lattice: weight 1 to an edge without 
an arrow, and weights u~ x / 2 , w^ 1 ^ 2 or u _1 / 2 iy -1 / 2 to a horizontal, vertical or oblique edge having 
an arrow. If needed, the edge weights can be included in the vertex weights, but we will not do 
that here. 

The statistical sum S(u, w) of our 6- vertex model is, of course, a sum of products of vertex and 
edge weights over all arrow configurations. The next lemma is the key statement. 
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Lemma 13.2. The statistical sum S(u,w) is a sum over the same sets of trajectories as the 
determinant I(u, w), and to each set corresponds the same summand up to, maybe, a minus sign. 
To be exact, the number of trajectories in the exponent of (—1) in Lemma 13.1 changes to the 
number of intersections (self-intersections included) of a given set of trajectories. 

Proof is evident from the statistical sum definition. 

Each closed path on the torus is homologically equivalent to a linear combination of two basis 
cycles a and b. The same is true for a set of paths (trajectories), regarded as a formal sum of 
them. Different sets may be homologically equivalent to a given cycle Za + mb but, as the following 
lemma shows, they have something in common. 

Lemma 13.3. For any set of trajectories on the torus homologically equivalent to a cycle Za + mb 
(a, b being basis cycles, l,m — integers), 

(number of intersections) — (number of trajectories) = Im — I — m(mod 2). (13.4) 

Proof may consist in the following simple consideration: 1) if the set consists of I trajectories 
going along a, and m ones going along b, (13.4) is obviously true, 2) under deformations of 
trajectories, the number of intersections changes only by even numbers, 3) with elimination of an 



intersection 



or inverse operation, the LHS of ( 13.4 ) may change also only by 



an even number. Starting from an arbitrary set and applying the transformations 2) and 3), one 
can arrive at a set of type 1), so the lemma is proved. 

Let the following products of edge weights correspond to the basis cycles: x — u^w^ 1 for a 
and y = u a2 w^ 2 for b. Denote 

s(x, y) = S(u, w), f(x, y) = I(u, w). (13.5) 



Theorem 13.1. The statistical sum of the inhomogeneous 6-vertex model on the kagome lattice 
defined in this section is invariant with respect to the evolution of the reduced 2 + 1-dimensional 
model (for all u,w) and is connected with the determinant I(u, w) ( 13.2\ ), whose vanishing defines 
the invariant curve of the model, by relations (in the notations of ( 13. S^) ) 



s(x, y) = 1/2 (-f(x, y) + f(-x, y) + f(x, -y) + /(- 



-y))-. 



f(x, y) = 1/2 (~s(x, y) + s(-x, y) + s(x, -y) + s(-x, -y)) . 



(13.6) 
(13.7) 



Proof. It follows from Lemmas 13.1, |13.2 , 13.3 that in the expansions of s(x,y) and f(x,y) in 
powers of x and y the coefficients near x l y m coincide if l m — l — m is even, and differ in their signs 
in the opposite case. This is exactly what the formulae ( 13.6] . 13.7 ) are about. 



Remark 13.1. From the physical viewpoint, the determinant I(u,w) can be regarded as a sta- 
tistical sum of the same model, but with other boundary conditions than for S(u, w). This can be 
explained as follows. Let us regard our torus as obtained by indentification of "opposite sides" of 
a plane domain that can be obtained by cutting the torus along the basis cycles a and h. We will 
consider those cycles as components of the plane domain boundary. Of course, we can determine 



the numbers I and m, for each set of trajectories from Lemmas 13.1, 13.2, 13. c, if we just know 



how the trajectories intersect with the boundary (and we may know nothing about the behavior of 
trajectories inside the domain). And it is exactly the numbers I and m that determine, according 
to the proof of Theorem 13.1, whether a sign must be changed of a given summand in a "usual" 
statistical sum S(u, w) in order to get I(u, w). 
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14. A finite-gap theory for the model on infinite lattice 



We considered in Section [13] an inhomogeneous 6-vertex model satisfying the "free fermion" con- 
dition, on the kagome lattice. The lattice itself was situated on a torus, thus, in particular, it was 
finite. The evolution was in fact "local" , i.e. a given weight in a moment r was influenced only by 
weights in a few neighboring points in the moment r — 1. However, the evolution was described 
"globally" , just as a particular case of the evolution from Section [?]. 

In this section, we will describe the evolution in local terms, the description being valid, be- 
cause of its localness, for an infinite in both spatial directions lattice as well. The evident solitonic 
character of the model stimulates one to study, in particular, the "finite-gap" quasiperiodic solu- 
tions on this infinite lattice. Recall however that the algebraic curve appearing in Section [l^ was 
associated with a "global" object — the statistical sum of the model regarded as a function of two 
parameters, u and w. In this section we develop another approach, namely, in the spirit of the 
usual theory of finite-gap solutions, we start from a given algebraic curve with some marked points 
in it and construct a solution out of those objects. 

It is natural to begin with a definition of evolution of matrices corresponding to vertices of a 
finite or infinite kagome lattice independent of any global objects. More precisely, we will assume, 
to be in accord with the evolution definition given in Section [7| that before a step of evolution 
there are given some products of such matrices (of size 2x2), while the matrices themselves appear 
and disappear within the step. Hence, at the beginning of a step, to each triangle of the form 



contained in the kagome lattice (and shaded in Fig. 3.7 situated on p. 38) a 3 x 3-matrix 

of complex numbers must correspond. The step begins with this matrix being factorized into a 
product of three matrices of the following form: 





h 




( a 2 





b 2 


C\ 


<h 


• 





i 


















d 2 




(14.1) 



Tins done, we will assume thai in our Irianjde the matrix ( Ul ^ ) corresponds to flic verlex 



the matrix ( ° 2 I corresponds to the vertex 
C2 d 2 ' 



, and the matrix 



03 

C3 



d 3 



corresponds to the vertex 



The kagome lattice can also be viewed as made up of triangles 



of the form 




(non-shaded in Fig. 3/7). Suppose that in such a triangle to the vertex 

c-2 d 2 



corresponds a matrix I ^} - 1 I , to the vertex 
1 ci d 1 



-a matrix 



-a matrix 



b 3 
d 3 



. Then let the matrix 



and to the vertex 




(14.2) 



correspond to the triangle as a whole. Finally, turn inside out all the triangles of the form 



in our lattice, converting them into triangles of the form 





(e.g. moving all oblique lines 



one lattice period right), putting the same matrix (14.2) in correspondence to each "converted" 
triangle as to the initial one. The description of a step of evolution is over. 

The same can be presented in a somewhat other way, starting from the triangular lattice of 
Fig. [3~3 (p. pq). By the beginning of an evolution step, to each circle — a vertex of that lattice — a 
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3 x 3-matrix must correspond that is to be factorized in a product ( 



14.1), which corresponds to 



a "decomposition" of the circle in a triangle 



with the triangles 



In these terms, the step of evolution ends 



being "packed" in new circles. 

Just as in Section [?], the evolution is defined up to "gauge transformations" . To explain what 
the gauge transformations look like now, note that each entry of the abovementioned 2 x 2-matrices 
naturally corresponds to a certain pair "incoming edge, outgoing edge" for the given vertex (Fig- 
will remind which edges are called incoming and outgoing) . If we mark those pairs 

' a\ bi 
d di 



ure 3.6 



on p. 



of edges with arrows, we will obtain the following pictures for entries of the matrix 



ai 



Cl 




h : 



di 




Now it is seen that to each edge a set of elementary gauge transformations naturally corresponds 
parameterized by the set of nonzero complex numbers. Namely, we can multiply by a nonzero 
constant the matrix entries in the vertex corresponding to the given edge as an outgoing one, and 
divide by the same constant the entries in the neighboring vertex corresponding to the given edge as 
an incoming one. Nonzero numbers can be put in correspondence to all edges, and the elementary 
gauge transformations can be performed for all simultaneously. The resulting transformation will 
be called simply a gauge transformation on the whole lattice. 

Lemma 14.1. 1 ) The result of the step of evolution described earlier in this section is determined 
to within a gauge transformation. 

2) The result of the evolution step is not changed if a gauge transformation is applied to the 
initial state. 

3) For a finite-dimensional system, the evolution coincides with that defined in Section [^j, and 
the gauge transformations considered here are such in the sense of Section |^ as well (but there are 
more gauge transformations in Section^ because they include also "non-local" , from the standpoint 
of model on the kagome lattice, transformations). 

Proof follows immediately from the definitions given here and in Section ^. 

Passing on to the "local" algebro-geometrical description of evolution, let us consider first the 
following abstract divisor evolution on an infinite in both spatial directions lattice. At the moment, 
we don't need to know the exact structure of those divisors or an algebraic manifold where they 
belong. Thus, let us temporarily understand by "divisors" just elements of some abelian group Q. 
Let six elements X>i, . . . , T>q be fixed in that group. 



Let a divisor correspond to each edge of the triangular lattice in Fig. 3.5 by the beginning of a 
step of evolution, with the following condition fulfilled: for each lattice vertex (a circle in Fig. |3.5| ) 
an element T> £ Q can be indicated such that the divisors corresponding to edges abutting on that 



vertex are as shown in Fig. 3.8 



Thus, T> depends linearly on the coordinates of a vertex, increasing by T>\ + T>2 — £*4 — £>6 when 
moving one lattice period to the right, and by T>^ + P4 — T>^ — T>2 when moving one lattice period 
up. 

A s tep of evolution begins with every circle being decomposed into a triangle such as depicted 
in Fig. 3J), where the numbers near an edge serve as a brief nota tion for a divisor corresponding 
to it, e.g. "12" denotes V — T>\ — T>2, and so on, compare Fig. \ZA\ Thus, the kagome lattice arises 

in it and "pack" 



from the triangular one, and then we take the triangles of the form 

them in circles — the vertices of the new triangular lattice. With this, the divisor evolution step is 
over. As the following lemma states, everything is ready for a next step. 
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Lemma 14.2. The divisors corresponding to edges abutting on each vertex of the triangular lattice 



obtained in the end of the previous paragraph are again of the form of Fig. 3.t , with D depending 
on the vertex. 

Proof follows from an easy direct calculation. 

Now let an algebraic curve T$ of genus go be given, and the abelian group Q let be the group 
of all divisors on r . Let T>i, . . . ,T> & be divisors on r , consisting each of one point (note in 
parentheses that, of course, one can consider a "vector" model as well, T>x, . . . ,V 6 consisting in 
that case each of several points). Let divisors V corresponding to vertices of triangular lattice 
(see Fig. [3~S| ) be of the degree go + 2. If all those algebro-geometrical objects are generic, the 
Riemann-Roch theorem shows that to each edge of triangular (Fig. |3.8| ) or kagome (Fig. |3.9D 
lattice corresponds a one- dimensional space of meromorphic functions / such that 

if)+V-V j -V k <0 ) (14.3) 

if a divisor V — T>j — T>k corresponds to the given edge. 

It remains to explain how the matrices are constructed corresponding to vertices of our two 
lattices. Fix a nonzero function / satisfying (14.3) for every j and k, and denote it as fjk- Then, 
say, the matrix -4i OC ai corresponding to a vertex in Fig. 3.8 is found from the relation 



A 



local 



fl2 \ / fi6 

/l3 J = I /56 
Jul V /52 



(14.4) 



that must hold in each point oftho curve r u - Similarly, the matrix ( Ql ^ ) (see (14. 1) and the 



graphical analog of that relation, Figure 3.9) is found from the relation 



ai h 
ci di 



he 



fi6 
/56 



(14.5) 



etc. Relations of type (14.4-14.5) determine the matrices correctly due to the fact that all mero- 
morphic functions in every such relation lie in a linear space of needed dimension. For example, in 



relation (14.5) divisors of all functions satisfy the condition 

Ujk) + V-Vo>0, 

hence the Riemann-Roch theorem, together with equalities 

deg£> = .9o + 2, deg£> 6 = 1, 

shows that the space of such functions is two-dimensional. 

Note that other choice of any function fjk, i.e. its multiplication by a nonzero constant, corre- 
sponds to an elementary gauge transformation of matrices. 

Lem ma 1 4.3. The described above divisor evolution generates, by means of formulae of type 
(14.4 ) and (14.5 ), the matrix evolution described in the beginning of this section. 



Proof. We must show that the operation of decomposition of a triangular lattice vertex in a 
triangle of the kagome lattice, and the inverse operation of "packing", when applied to divisors, 
generate the similar operations on matrices. E.g., if a matrix ^4i OC ai is determined by the relation 
(14.4), and we have to factorize that matrix in a product (14.1), then the values a\, . . . , d^ obtained 
from (14.5) and two similar relations must give the solution to the problem. It is worth while to 
write out explicitely those two relations similar to (14.5): 



a 2 b 2 
C2 d 2 



fl2 
/23 



/26 
/52 



0,3 b 3 
C3 d 3 



fl3 
/34 



/36 
/23 



(14.6) 
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If we now let the product (14.1) act on the column from the LHS of (14.4), and use (14.5) and 



(14.6), we will get exactly the column in the RHS of (14.4), which means that the factorization of 
a matrix is done with the help of divisors in our local case with equal success as in the global case 
of Section fTo], cf. formulae (10.4-10.6). Certainly, the similar considerations are valid also for the 
"packing" operation. The lemma is proved. 

Thus, we have shown how to construct solutions to the problem of matrix evolution on the 
infinite lattice in the spirit of standard finite-gap integration theory In Section ^ we will connect 
this approach with the "global" approach of Sections |9| and [lOl Explicit expression for the solution 
in multidimensional theta functions is presented in Section hq. 



15. Connection between "local" and "global" curves in the 
periodic case 



We continue to consider divisors from Section 14 corresponding to edges of the infinite triangular 
lattice. Let there exist integers £i and iji such that the divisor equivalence holds in the curve Tq 



£i(X>! + V 2 - P 4 - X> 6 ) + »7i(X>3 + ^4 - V 5 



0. 



(15.1) 



(compare with the text between lemmas 14.1 and 14.2). As each of the divisors T>\, . . . , Vq consists 
now of just one point, it is convenient to introduce for those points notations Pi, . . . , Pq. Note 
also that V£s in ( |15.1 ) arc now other divisors than the "global" ones denoted by the same letters 
in Sections H-IO. Equivalence (15.1) means that there exists a function gi on T that has a zero 
of multiplicity £i in the point Pi , a zero of multiplicity £i — r/i in the point P 2 , . . . , a pole of 
multiplicity £i in the point Pg. Thus, meromorphic functions /, out of which we will construct 



matrices according to relations of type (14.4), can be put in correspondence to lattice edges in 
such a way that the multiplying of a function / by g\ will c orrespond to a lattice translation by 
the vector (£1,771). As for the matrices like -4i OC ai from ( 14.4 ), they will obviously be 'periodic with 
period (£i,?7i). 

Lemma 15.1. Let there exist two linearly independent vectors with inte ger e ntries (£1,771) an d 
(£2,772) such that two following divisor equivalences in the curve Tq hold: (15.1) and a similar one 



£2 (X>i + V 2 - V 4 - X> 6 ) + 77 2 (p a +V 4 -V 5 



0. 



(15.2) 



Then matrices corresponding to lattice edges are, in a proper gauge, doubly periodic in coordinates, 
with periods (£1,771) and (£2,772). 

Proof is obtained by adding to the above consideration for (£1,771) a similar consideration for 
(£2,772)- 



Thus, conditions (15.1) and (15.2) are sufficient for the model on infinite lattice to become, in 
essence, a model on a torus introduced in Section^. Our next task is to learn how to pass on from 
the curve To to the "global" curve T, and by that calculate the statistical sum of the model on 
the torus starting from a "small" curve Tq. To be concrete, consider the problem of constructing 

' X' 



vectors 



such that 









( uX 






- 


uwY 








\ wZ 



(15.3) 



the components of vectors X, Y, and Z corresponding to horizontal, oblique, and vertical edges 
respectively; A being a global operator, as in Section [13[ and u and w being some so far unknown 
values. 
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It turns out that the solution of (15.3) is obtained if u and w are introduced as multivalued 
functions on Tq by the formulae 

u^w^=g 2 , (i5 ' 4j 

where the function g 2 is constructed out of ( 15. 2| ) in the same way as g\ out of ( 15.1 ). The vectors 
X, Y and Z are built now as follows. Choose some vertex to be the origin of integer- valued 
coordinates (£, rf). For each of the vectors X, Y and Z, one componen t — t he one corresponding 
to the incoming edge abutting on the origin of coordinates (recall Fig. 3.6 on p. pq ) — is defined 
as the value of meromorphic function / corresponding to this edge as in Section [bijin some point 
of r . Next, we take some values u and w for this point according to (15.4), and for an incoming 
edge at a point with coordinates (£, 77) we take as a component of the vector X, Y or Z the value 
of the corresponding function / multiplied by u~£w~ v . As a result, as was required, the vector 
components are not changed under a translation by the periods (£1, J?i) and (£2, V2) (we assume of 



course that a "quasiperiodic" gauge was chosen for functions /, as in Lemma 15.1 



Lemma 15.2. The vectors X, Y and Z constructed in the previous paragraph satisfy indeed the 



equation (15.3) for any point of the curve T and any pair (u,w) lying above it (i.e. satisfying 



Proof. Indeed, it is easy to see that (15.2) is equivalent to the totality of relations of type (14.4) 
for all lattice vertices, because the factors u, uw and w in ( |15.3[ ) are canceled by factors arising 
because of multiplying the functions / by u~£w~ n . 

Having constructed the sol ution s of equation ( 15. 3| ), we arrive at some algebraic dependence 
between u and w arising from ( 15.4 ) and an algebraic dependence between the meromorphic func- 
tions 51 and g 2 on the curve Tq. We expect that it will be, maybe under some additional conditions, 
the same dependence as given by the "global" equations (9.1, 9.2), because solutions of ( 15.4 ) do 
satisfy, of course, the system (9.1, 9.2). However, we don't know at the moment whether our local 
construction does not lead to some singular operators A, e.g. those having an invariant curve that 
contains several components. In such case, one can conceive a situation where ( 15.4 ) together with 
the dependence between g\ and g 2 yields only a part of the curve T, while all T must be known, 
of course, for calculating the statistical sum. Besides, we must attach the exact meaning to the 
words "algebraic dependence between gi and g 2 " . 

It turns out that there exist simple conditions that guarantee the irreducibility of T. They are 
as follows: 

(i) irreducibility of To, 

(ii) the parallelogram built on vectors (£i,?7i) and (£2,772) must be the minimal parallelogram of 
periods (i.e. the parallelogram of minimal area among the parallelograms built on two linearly 
independent vectors of periods), and 

(iii) at least one of the pairs (£1,^2), (^1,^2), (£1 — ^1, £2 — V2) consists of relatively prime numbers — 
the condition without which the proof of the following lemma fails. 

Lemma 15.3. Consider the mapping <f> : Tq — > CP% x CP\ given by the formula 

z ^ (gi(z),92(z)). 

Under the above conditions (i) and (Hi), 4> JS a birational isomorphism of the curve Tq on its 
image <j)(T ). 

Proof. Compose the following table of zero and pole multiplicities for functions g\ and 172 in 
points Pi, . . . , Pq. 





Pi 


P2 


p 3 


Pi 


p 5 


p e 


91 




61 - m 




-ft + m 






92 


6 


6 - V2 


V2 


-6 + m 


-V2 


-6 
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This table is made up according to formulae ( |15.1| , |15.2| ). An integer in a table cell is the zero 
multiplicity of a function in a point if it is positive, or minus pole multiplicity if it is negative. 
Functions g\ and gi have no other zeros or poles. 

Suppose that the mapping is a g-sheeted covering. Let us replace 0(To) by its nonspecial 
model T ns which exists according to §7A of the manual [|o). Somewhat freely, we will denote the 
covering T$ —* T ns by the same letter 4>. 

Let us prove that a point Pj in which zero or pole multiplicities of functions 171 and gi are 
relatively prime (such a point exists according to the condition (hi)) cannot be a branch point of 
covering cj). Indeed, due to the indicated relative primality, a local parameter in the point Pj can 
be chosen as a product of some degrees of g\ and gi, i-e. is uniquely determined by an underlying 
point from a neighborhood of (j)(Pj), which cannot happen in a branch point. 

The point Pj from the previous paragraph also cannot have the same image as another point 
P G To, i.e. 4>(Pj) 7^ 4>{P) when P 5^ Pj. To show this, we note that P can only be one of the 
points Pi, ... , Pq, because only there zeros and poles of g\ and gi are situated. However, by virtue 
of the linear independence of periods (£1, 771) and (£2, 772), and the above table, the multiplicities of 
zeros or poles of gi and gi in points 4>(Pj) and <p(P) cannot coincide (even if P is a branch point). 
Hence 4>(Pj) ^ <j>{P). 

Combining the results of two preceding paragraphs, we see that Pj is the only point lying above 
4>{Pj), and not a branch point. Hence, the number of sheets q = 1, and the lemma is proved. 

We continue to consider the algebraic dependence between functions g\ and gi. In view of the 
irreducibility of Fo, this dependence can be expressed as 

Po(j9i,to) = 0, (15.5) 
Pq being an irreducible polynomial in two variables over the field of complex numbers. By 



Lemma 



15.3 



the curve given by equation ( 15.5 ) is birationally isomorphic to T. 



o- 



Theorem 15.1. Assume conditions (i), (Hi) on p. 46Connection between "local" and "global" 



curves in the periodic caselemma.l5..2. Then, equation (15.5), which describes the algebraic de- 
pendence between functions gi and g% realizing the equivalences (15.1, 15.2) of divisors on T , after 
substituting in it the expressions ( 15. 4 ) of gi and gi through u and w gives exactly the invariant 



curve of global operator A on a torus with periods (£1, 771) and (£2, 772) • 

Proof. Note that the equation of curve T has as well the form 

P(gi,g 2 ) = V(u^w m ,u^w n2 ) =0, (15.6) 

where gi and gi correspond to going around two basis cycles on the torus, see Section |l3|. We will 
prove that equations ( 15.5Q and fll5.6|) define the same curve. 



According to the table on p. [4q, functions gi and gi have li and I2 poles respectively in Tq, 
where 

l j =2-ma X {\Z j l\r) j l\Z j -r lj \}, j = l,2. (15.7) 

Hence any generic value is taken by a variable gj in lj points of the curve, so that lj values of the 
other variable correspond to it. This situation can be described by equation ( 15. 5| ) only if Vq h 
degree li in gi and degree li in gi 



as 



Passing on to equation ( 15.6 ), note that it is natural to suppose that the degrees of gi and 
gi (but not u and w) in it can be negative, if the expansion of the cycle associated with a given 
trajectory (see Section |l3|, three paragraphs starting from the one containing formula ( 13. 2\) ), in 



basis cycles has one or two negative coefficients. It is not hard to check that the difference between 
the maximal and minimal degrees of, say, gi docs not exceed the number of lines in the triangular 



lattice of Fig. 3.5 that intersect with a vector with coordinates (£1,771) (it is convenient that we 
imagine it as situated "generically" , e.g. having a point with irrational coordinates as its origin). 
Indeed, the difference between the greatest and least intersection numbers of any trajectory going 
along the edges on the torus always to the right, upwards, or to the north-east, and the basis cycle 
corresponding to the vector (£1, 771), cannot exceed, of course, the mentioned number of lines. This 
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number of lines is, however, nothing else but li from formula (15.7), because the vector (£1,771), 
evidently, intersects with |£i| vertical, \r]i\ horizontal, and |£i — 771 1 oblique lines. 

Thus, we see that, having multiplied V(g\,g-i) (if needed) by its common denominator, we get 
a polyno mial o f the same degr ees in variables 171.2 as Vo(gi, 92)- Hence it is clear that the 



equation ( 15.5 ) together with ( 15.4 ) yields indeed the whole "global" curve T. The theorem is 
proved. 

The following lemma describing more precisely the structure of T will be the last in this section. 
It is here that we will use the condition (ii), not needed for us before. 

Lemma 15.4. Consi der a covering T' of the curve Tq defined by multivalued functions u and 



w according to ( 15.4 ). Let the conditions (i) and (ii) on p. 46Connection between "local" and 
"global" curves in the periodic caselemma.l5..2 hold (while we will not use condition (Hi) that 
guarantees that T' and T are isomorphic). Then the curve T' is irreducible. 



Proof. The system (15.4) defines an extension of the field of meromorphic functions on Tq. 
This extension is normal, because there exist automorphisms of T' mapping its any sheet into any 
other, defined by an evident formula 

(u, w) — » (u}iu,u) 2 w), (15.8) 
oj\ and UJ2 being some root of unity. There are |£i?72 — £,2Vi\ such automorphisms, and if they all 



enter in the Galois group, then this group acts transitively on solutions of the system (15.4), and 



both this system and, hence, the curve V are irreducible (see a chapter on Galois theory in the book 



|79[ or ]84|]). In this case, besides, the extension defined by system (15.4) is | ^ 1 772 — £2771 |-dimensional 
over the function field on Tq. 

If, otherwise, the curve V is reducible, then this extension is less then |£i?72 — £2771 |-dimensional 
over the function field on Tq. Consider all possible products u^w v of integer degrees of u and w 
for £ and 77 lying within one parallelogram of periods (it is implied that every such parallelogram 
contains |£i?72 — £2771! integer points, and that one can cover all the plane with such parallelograms, 
with periods (£1,771) and (£2,772), while other details of their disposition are not important). The 
products vfiw 11 are multiplied under the action of Galois group elements by some characters of that 
group, and as there are |£i772 —£2771! products, and fewer characters, then within one parallelogram 
of periods there are products to which correspond the same character. This makes clear that there 
exists a product go — u^"w n " invariant under the Galois group for which (£0,770) is not an integer 
linear combination of periods (£1,771) and (£2,772). The existence of such function go on Tq means, 
however, that (£0,770) is a period, which contradicts to the minimality of parallelogram built on 
(£1,771) and (£2,772). Thus, an assumption of reducibility of T' has led to a contradiction. The 
lemma is proved. 

Let us sum up the results. In this section, a situation was described when matrices in vertices 



of the triangular lattice constructed by method of Section 14 depend doubly periodically on the 
coordinates, so that a model on an infinite lattice becomes a model on a torus. For such a model, 
a statistical sum exists that is a polynomial in u and w closely connected with the equation of 
invariant curve of a "global" matrix A composed of matrices situated in lattice vertices. To deal 
with this situation, we presented, under some restrictions, a principle of constructing a "global" 
invariant curve (usually, of a very high genus) out of a curve Tq (that may be, e.g., elliptic). The 
mentioned restrictions guarantee the irreducibility of algebraic curves whose components otherwise 
might be lost or confused. A constructive realization of the mentioned principle will be presented 
m Section [l7[ 

Note that the statement about obtaining the whole T was deduced without using condition (ii) 



and Lemma 15.4. We used only the irreducibility of V as a polynomial in g± and 52- Condition (ii), 
guaranteeing the irreducibility of V as a polynomial in u and w, will not be needed for the ex plicit 
calculation of statistical sum in Section [Tt]. It appears however that a not complicated Lemma 15.4 



gives a useful information on the curve T. Later we will use this lemma to calculate the number 



of components of some reducible curves T (Remark 17.2 ) 
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16. 



Expression of a finite-gap solution in theta functions 
and a multilinear form of the equation 



In this section, we will express the meromorphic functions fjk entering in formulae ( 14.4 - [l4~6|) and 
defined on the curve Tq, through a theta function 8(z) defined on the Jacobian Jac(To) of that 
curve. Recall (see §1 in Chapter II of the manual |8l|) that 



0(z) = 0(z, fi) = E exp(7rm T fin + 27rm T z), 

n G ZBo 



(16.1) 



where z is a complex column vector of height go equal to the genus of curve To; the components 
of vector n of height go run through all integer values; and Q is the matrix of periods of To- We 
will obtain, from the expressions for fjk, very simple expressions for matrix elements of matrices 
entering in LHS's of ( 14. 5| , 14.6 ). This, in its turn, will lead us to constructing of some analogs 
of Hirota's r-function satisfying a six-linear homogeneous equation. As is known, recently there 
appear more and more integrable equations for which one cannot construct a bilinear Hirota 
representation. It is replaced, according to papers |87]], by its multilinear analog. 

Pass on to a detailed account. As was explained in Section [l4|, meromorphic functions corre- 
spond to edges of the kagome lattice, and out of those functions matrices can be found correspond- 
ing to vertices of type , or S , according to the formulae ( |14.5| , |14.6 ). We will 

write out a general expression for such functions through the theta function ( |l6.l[ ). Choose one 
of the lattice edges as "initial" and divide the meromorphic functions on all edges by the function 
corresponding to the initial edge. Hence, we change all divisors satisfying relations of type 



(/on a given edge) ~f~ ^on a given edge _\ 



(16.2) 



by a constant divisor. 

Now to the initial edge a function satisfies identically equal to unity, and also an effective divisor 
"Dini, i.e. a formal sum of go points in r . According to §3 in Chapter II of the book |^5[ , we can 
choose a vector Zq G C 90 and a point Pq G To in such a way that the multivalued function 



<p(P) = 9( 2o + J < 



(16.3) 



of a point P G To will have exactly Z?i n i as its zero divisor. Recall that uj is a vector composed of 
holomorphic differentials on IV 

As for the rest of the edges, to each of them a divisor of type 



V, 



on a given edge 



k k 

E^-E^ 

3 = 1 j=l 



(16.4) 



corresponds, where k is some number (increasing as we are moving off the initial edge), and each 
of the points Rj and Qj coincides with one of the points Pi, ... , P§. Correctness of (16.4) follows 
immediately from comparing divisors on neighboring edges, see e.g. Fig. [3^ (from where it is seen 
also what the points Rj and Qj exactly are; we will be concerned with this somewhat later). 
The explicit expression for a function satisfying (16.2) with a divisor (16.4) is as follows (nec- 



essary explanations are given just below the formula): 

e(z + /w-E 

Jon a given cdgc(^) 

eUo 



Qj 



3=1 



Qj 



Po 



k 

n< 



(16.5) 
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P2 


+ Pe- 


-P3-P4 






ls 


k 



P2 + P6-P1- P3 



Figure 3.10. The "initial" edge and divisors /J-Rj — J]] Qj on the neighboring edges 

3=1 3=1 



Here e £ C 90 is any such vector that 



0(e) = 0, 



(16.6) 



while the function 



.v 

E e (x,y) = 9(e + J u) 



of two points x,y GT, called principal form, does not identically vanish. It remains to indicate the 
integration paths in (16.5), more precisely, their homological classes. Connect all points Rj and 
Qji 1 < j < k, and also P, with Pq by arbitrary paths. This done, we will assume that each path 
QjRj consists of two parts QjPo and PoRj, and, similarly, that paths RjP and QjP also "pass 
through the point Pq" • 



One can verify by standard means that (16.5) defines a single- valued function on the curve To, 
so that the choice of path PqP in fact plays no role. As for the equality ( 16. 2| ), it follows from 
the fact that all "superfluous" zeros of principal forms in the enumerator and denominator of the 
right-hand fraction in (16.1) cancel each other p5[ . 

Now pass on to calculating matrix elements corresponding to transitions between neighboring 



edges of kagome lattice, by formulas like (14.5, 14.6), indicating en passant concretely points Rj 
and Qj in formulas ( |16.4| , |16.5 ). To begin, consider a situation depicted in Fig. 3.1C. This figure 
is, in essence, a modification of Fig. 3.9. The thick edge is the initial one. Near other edges, the 

k 



k 

divisors \_. Pj ~ 

3=1 j=i 
are situated in points R 



Qj corresponding to them are written out (recall that the poles of functions 



wh ile the ze ros — in points Qj). The system (14.5), after substituting the 
functions according to ( 16. 5|) and Fig. 3.10 and multiplying by the common denominator, acquires 
the following form: 

/ P P \ 

0(z o + fu) 0(e- 



ai 

Cl 



h 
di 



P-2 



9(z 



Po 



P:< 



P-, 
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e(z 



p p 2 

u> — I U}\ 



Pa 
Pi 



e(z 



V 



lj — 



6 e 



die 



Pi 



(16.7) 



Ps 



} 



Po Ps 

The equality (16.7) must hold for any point P G IV Recalling ( 16. 6| ), we see that a\ and C\ 
are easily found if we set P = P 3 , while b\ and d\ — if we set P = P2 . It is convenient to introduce 
notations 

p p 3 



z{P) 



1 < j < 6. 



Po 



We get 



ai = 



6>(z + z 3 - z 2 + z 4 ) 0(e + z 3 - z 4 ) 



(z 
(zo 



- z 3 ) ! 

z 4 )( 



le- 
fe- 



- z 3 - z 2 ) 
z 2 - z 4 ) 



ci 



0(z o +z 3 )0(e- 
6»(z + z 3 - z 2 + z 5 



z 2 - z 3 ) ' 
9(e + z 3 - z 5 ) 



(z +z 3 )6»(e- 
'(z +z 5 )6»(e- 



Z3 - z 2 ) 
- z 2 - z 5 ) 



(16.8) 
(16.9) 
(16.10) 
(16.11) 



Matrix elements a 2 , . 



6»(z + z 3 ) 0(e + z 2 - z 3 ) 
, (i 3 are found similarly. We will write out for them only systems of 



type (16.7), because explicit expressions like (16.8-16.11) are obtained from them in an obvious 
way by using ( 16. 6| ). So, 



"3 
C3 



02 
C2 



b 3 
d 3 



62 
d 2 



( 6{z + z{P) - z 6 + zi) 6»(e + z(P) - zi) 
I 0(z o + z(P) - z 6 + z 3 ) 0(e + z(P) - z 3 ) 



0(z o +z(P))0(e + z(P)-z 6 ) 
9(z + z(P) - z 6 + z 5 ) 0(e + z(P) - z 5 ) 

zo + z(P) - z 2 - z 6 + z 3 + z 4 ) 0(e + z(P) - zi) 
>o + z(P) - z 2 - z 6 + z 3 + z 4 ) 0(e + z(P) - z 4 ) 



/0(z o + z(P) -z 2 + z 3 ) 
\e{z + z(P) -z 6 + z 3 ) 



le- 
(e- 



z(P) 
z(P) 



Z6; 

z 2 ) 



(16.12) 



(16.13) 



No we can extend Figure 3.10 in any direction and calculate matrix elements in some more 



vertices. Fortunately, not complicated calculations like those done above show that the formulae 



(16.7-16.15) remain valid under translations by lattice periods, if a change 



zo 



zo 



£(P 4 + P 6 - Pi - P 2 ) + V (P 2 + P 5 -P 3 - Pi), (16.14) 
is done in them, where the integers £ and r\ show show how many lattice periods horizontally and 



vertically we have moved off the "initial" triangle in Fig. 3.1C. 

Perfectly similar considerations, with using the divisor evolution from Section |l4|, are valid for 
the dependence of matrix elements on the time. The result consists, of course, in adding in (16.14) 



a linear dependence of z on the time r. We will formulate this as follows: if divisors on edges 

1 1 

then 



around two vertices of the same type (for example, of type —/— ) differ by Rj — Q 



j ■ 



the corresponding values z differ 



1 Vj 
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Let us end this section by some simple remarks about an analog of Hirota's bilinear representa- 
tion for our dynamical system. As was indicated in the beginning of this section, only a six-linear 

representation could be found, which is obtained if we write the entries of matrices like ' ' 1 1 

in the form 

a\ = — , h = — , ci = — , di 



ti ' 1 ~ h ' 

etc. and then multiply every equality of the form 



7i 



Si 



ci di 
(16.15) 




(16.16) 



(where the "shaded" matrices belong to the time less by a time unit than the non-shaded ones) by 
the common denominator t'^t'^t^ts. The reasonableness of the obtained equation is substantiated 
by the following. Recall the formulae ( 16.8 - 16.11 ), where only Zq depends on coordinates and time, 
and set ti = 8(z +z 3 ). Similarly, choose the rest oft , equal to theta functions of z +const situated 
in the denominators of formulae of type ( 16.8 - 16.11 ). We get for each of the values a j , j3j , jj , Sj,tj 
that it is in the "finite-gap" case a theta function multiplied by a harmless constant, which is a 
distinguishing feature of a Hirota r-function. 



17. Invariant curve equation in the thermodynamical limit 

As soon as we began in Section [li] to consider an inhomogeneous 6- vertex model on the infinite 
lattice whose Boltzmann weights are constructed out of an algebraic curve Tq (of not very high 
genus) and some divisors in it and depend on coordinates quasiperiodically, we naturally come to 
the problem of existence and calculation of the thermodynamic limit of specific free energy, i.e. the 
limit of the logarithm of statistical sum divided by lattice area. In this work, some steps in this 
direction are propose d. Ju s t now we will continue to consider the doubly periodic case of Section [L5[ 



assuming conditions ( 15.1 , 15.2 ) of linear dependence with integer coefficients for divisors Pi — Pq, 
Pi — P4 and P 3 — P 5 . The vectors of periods (£1, 771) and (£2, ^2) can be, however, as big as desired. 
We will see that there exists an elegant integral representation for the real part of the limit of 



the ratio of the logarithm of determinant ( 13. 2| ) to the area of the lattice "covering" the lattice on 



the torus with periods (£i,?7i) and ( £2, rh ), with both dimensions of the covering lattice tending 



to infinity. It follows from Theorem 13.1 on the "physical level of rigor" that this limit coincides 



with the specific free energy (see also Remark 13.1). We will present some confirmation for this 



in a simple particular case in the concluding part of this section, while now returning to the 
mathematical level of rigor. 

Let meromorphic functions g{z) and h(z) be given on a smooth irreducible algebraic curve Tq. 
Then those functions obey an algebraic dependence 

V(g,h) = 0, (17.1) 

V being an irreducible polynomial. Standard methods of the theory of functions of a complex 
variable provide in a not complicated way an integral representation for the absolute value of 
V(g,h). 

Introduce the following notations: k and I will be the numbers of poles, with regard to their 
multiplicities, of functions g(z) and h(z) respectively. By gj(ho), where 1 < j < I, we will denote 
the values of function g in those points where h(z) = ho (we take those points in any order and 
also with regard to their multiplicities). Similarly, hj(go), with 1 < j < k, will be the values of h 
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in points where g(z) = go- The following integral representation (a version of the Cauchy integral 
formula) takes place: 



2 J> 

3=1 



9o-gj(ho) 



go -5j(oo) 



= 2 £ln 

3=1 



h - hj(go) 



ho — hj(oo) 



1 

2tH 



dg(z) A dh(z) 



l{z) - 9o)(h(z) 



(17.2) 



Let us prove, for the sake of completeness, the formula ( 17.2 ). Consider a real- valued function 

<p(z) = 2ln\g(z) - g \ 

dehned on To, except, of course, the points where g{z) = go or oo. As the holomorphic part of the 
function ip differential coincides with dg(z)/ (g(z) — go), the integral in (17.2) can be rewritten as 



1 

2tH 



dtp(z) A dh(z) 
h(z)-h 



(17.3) 



By Stokes theorem, the integral ( 17.3 ) equals 

97T7 ^ 



over singular 
points 



ip{z) 



dh{z) 



h(z) 



where the sum is taken over all singular points of the integrand in (17.3), while the integrals under 
summation sign — along infinitesimal contours around those points in the anticlockwise direction. 
As the singularities of function ip have the character of modulus of local parameter logarithm, 
the integrals along infinitesimal contours around them vanish unless the singularities of ip coincide 
with the points where h{z) = ho or oo. The integrals around those points, divided by — 27ri, equal 
2 In \ gj(ho) — go\ and 2 In \gj (oo) — go\ respectively. Hence 



2$> 

3=1 



go - 9j{ho) 



go - g 3 {oo) 



1 

2ttI 



dg(z) A dh(z) 



(g(z) - 9o)(h(z) 



The remaining equality in (17.2) follows from the obvious equalness in status of functions g and 
h. 



To link the equalities (17.2) to the polynomial V(g,h) from formula (17.1), note that 



V(g,h)V(g',h') 



lim , , 

3 '^oo V(g', h)V{g, hf) g - 9j (oo) 

h'-*cr 



n 



9~9j(h) 



n 



h-hj(g) 



(17.4) 



h — h-i(oo) 

Equalities ( |17.4| ) are proved by directly comparing the zeros and poles of LHS and RHS's in variable 
g or h. Let V(g, h) be normalized so that the coefficient at the term of highest degrees in both 
variables, i.e. g l h k , equals unity. We have 



v(g',h) 
(g'Y 



no-^(oo)), 

3 = 1 



(17.5) 
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(h') k 



n(3~5j(oo))- 



Formulae (17.2) and (17.4-17.6) together yield 



(17.6) 



\V{goM)Y 



exp 



2m 



dg(z) A dh(z) 



(g(z) - 9a)(h(z) - h ) 



Y[( h ° ~ h o(°°)) ■ IIC9 -ffp(°°)) 



P =i 



(17.7) 



We will see later that, in order to pass to a "thermodynamical limit" in formula ( |17.7[ ), we must 
consider the following situation. Let functions g(z) and h(z) on T be represented in the form 



h(z) 



(17.; 



where £ and 77 are integers, while u and w are functions on some covering of To (this can always 
be done, of course). Let a and b be positive integers divisible by £ and 77 respectively. Introduce 
functions 



G(z) 



a/f 



H{z) 



w 



h(z 



Consider a limit 



where 



1 



b— »OG 



dG(«) A dH(z) 



^yy (g(z)-g )(h(z)-h )' 
1 



G 



(17.9) 
(17.10) 

(17.11) 



with some Uq and Wo n °t depending on a and b. Changing the order of passage to limit and 
integration (the legitimacy of which operation is verified by standard means), we find that the 
limit ( 17.10 ) equals 

J J dlnu(z) Adlnw(z), (17.12) 
A 

A being that part of To where 

|u| > K|, H > Kl- (17.13) 

Let us now link the above considerations to Section [l^ in the following way. Let us take 
functions gi(z) and g%(z) from Section ^5|as g{z) and h(z) respectively. Assume for simplicity that 
111 formulae ( jl5.l[ ||J) 

Vi =6 = 0, (17.14) 



and denote 



v = m- 



Thus, the torus becomes "rectangular", with periods (£,0) and (0, 77) , while formulae (15.4) turn 
into (17.8). Assume also conditions (i-iii) on p. 46Connection between "local" and "global" curves 
in the periodic caselemma.l5..2, observing that the last of them will mean just that £ and 77 are 
relatively prime. 



It follows from the table on p. [16| and formula (17.14) that g\{z) has zeros of multiplicity £ in 
points Pi and P2, and poles of multiplicity £ in P4 and Pe, while 92(2) has zeros of multiplicity rj in 
points P3 and P4, and poles of multiplicity 77 in P2 and P5. Hence, it is clear that functions u and 
w on a ^77-sheeted nonramified covering of Tq possess the following zeros and poles of multiplicity 
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zeros 


poles 


u 


Pi, Pi 


Pi, Pe 


w 


P 3 , Pi 


P2, P 5 



Remark 17.1. Functions u and w, thus, can be expressed in a simple way through principal forms 
on Tq Ka/, but we will not need those expressions here. 



As will be clear soon, the limit 



lim — In 

f-* 00 ab 

b — >oo 



V(G ,H ) 



is important for us. Here the polynomial V(G, H) defines the algebraic dependence between func- 
tions ( 17.9Q which can be written as 



V(G,H) = jil[v( 9 exp 2 ^kexp^ 



p=l q=l 



0. 



(17.15) 



Substituting G and H in place of g and h in (17.7), taking into account the equality of expressions 
(17.10) and (17.12), and knowing the poles of functions G(z) and H(z) (see ( 17.S| ) and the table 
before Remark 17.1), we find: 



lim — In 

ab 

b — s-oo 



V(G ,H 



^— // d lnu(z) A dlnw(z) + In |u | + In |wo| + 



+ ln(max{K|, |«(P 5 )|}) + ln(max{K|, |w(P 4 )|}) , (17.16) 

where domain A C Tq is, as before, defined by inequalities ( |17.13 ), while uo and wo, of course, 
obey ( |17.11[ ). 



Pass now on to explaining why the polynomial (17.15) arises when we take a lattice on the 
torus a x b that covers torus £ x r\. 

Theorem 17.1. Let the invariant curve corresponding to the model on "rectangular" torus of 
sizes £ x rj be given by equality 

r(g,h) = V(v,t,w r >) =0. (17.17) 

Let numbers a and b be divisible by £ and n respectively. Then the invariant curve corresponding 
to model on the lattice on torus a x b, with matrix elements in the vertices "raised" from the 
lattice £ x 77 according to the natural covering, is given by equation ( 17. 1^ , where (17.8) must be 
substituted. 



Remark 17.2. Thus, the invariant curve in this case consists of copies of the covering of T 
dehncd by multifunction u and w ( |l7.^ j. Recalling conditions (i—ii) on p . 46C onnection between 
"local" and "global" curves in the periodic caselemma.l5..2 and Lemma 15.4 , we see that under 



these conditions the invariant curve contains exactly e components. 



Proof of Theorem 17. \ Recall that an invariant curve is birationally isomorphic to a set of 

such pairs (it, w) for which there exists a vector X = | Y ] such that 

Z 




I uX 

uwY 
\ wZ 



(17.18) 
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(see Definition 9.1 and formula (9.3) after it). If we deal with the torus a x b, then each of column 



vectors X, Y, and Z has ab components. Recall that when considering equalities of type (17.18), we 



pass from the kagome lattice to the triangular one, "packing" triangles in circles, as was explained 
in Sectio n |l3|. L et (x,y) be integer-valued coordinates of a given vertex of triangular lattice. The 
equality ( 17.18 ) remains valid under the transformation 



X 




exp(2iricix / a)X 
exp(2iricix / a + 2irc2y /b)Y 
exp(27ric2j/ /b)Z 



(17.19) 



(17.20) 



u — > exp(— 2irici/a)u, 
w — » exp(— 2nic2/b)w 7 

where c\ and ci are any integers. If matrix elements of A were raised from the lattice £ x -q, 
then among the vectors X there are periodic ones w ith periods (£,0) and (0, 77), and pairs (u,w) 
corresponding to such vectors form the curve ( 17.17 ). Clearly, all vectors X must be periodic with 
periods (a, 0) and ( 0, b), and they all can be obtained from vectors with periods (£,0) and (0, 77) 
by transformatio n (17.19] ) with proper c% and ci- It is not difficult to see that the corresponding 



transformations ( |l7.20[ ) of pairs (u, w) provide exactly || isomorphic components of the invariant 
curve according to ( |17.15[ ). The theorem is proved. 

Note that "integral of motion" I{u,w) ( |13.2| ) from Section [l^ differs from polynomial V(G,H) 
by a factor whose absolute value is \uw\~' Iab (see (17.15) and a remark about nomalization of 



V(g,h) before formula (17.5)). Hence we will get 



1 



1™, ~T ln l J Oo,W )| 
a->oo a0 
b — >oc 



if we subtract (21n|ito| +21n|w |) from the RHS of ( 17.16| ). Thus, the desired thermodynamic 
limit has been calculated. 



Remark 17.3. Logarithmic differentials of multifunctions u(z) and w(z) on the curve T , that 
occupy a key position in formula (17.16), astonishingly belong to the same type as differentials 
playing a key role in the algebro-geometrical string theory |g§[ |9(]/. Namely, they are differentials 
of the third kind (i.e. possessing only simple poles) on Tq, with all their periods pure imaginary 
(because absolute values of u(z) and w(z) are single-valued). The algebraic tech niques use d in this 
work forced us, however, to impose the integer linear dependence conditions (15.1, 15.2) on the 
divisors. 

Besides, formula ( \17.16\ ) reminds of string theory because the double integral in it can be 
interpreted as the area of domain A in Kahler metric corresponding to imaginary part of the 2- 
form dlnw(z) A dlnw(z) (changing the form to its imaginary part does not change the integral 
in (17.lt), because both sides in that formula are real). As is known, the "Nambu-Hara-Goto" 
action in string theory is proportional to the "world sheet" area \91 



In the remaining part of this section we will show how our method works in case of a homoge- 
neous model, choosing a curve of genus as Tq: 



CP 1 . 



(17.21) 



Clearly, conditions (15.1, 15.2) of divisor equivalence in this case always satisfied, and any integer 
vector is a period (in the sense of Lemma 15.1). 

To simplify the task still more, we will consider a model on square lattice. We will obtain it by 
"removing oblique lines" from the triangular lattice. For this purpose, we will set 



(17.22) 
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and see how this will affect the form of matrix -4i OC ai corresponding to a vertex of Fig. 3.8 type 
and determined from the relation (14.4). The latter is worth being written out here once more: 








)- 


/56 







(17.23) 



(recall that divisors T>%, . . . ,Dq on Fig. |3.8| consist each of one point Pi, . . . , Pq, and each mero- 
morphic function fjk equals zero in points Pj and Pk). 



Lemma 17.1. Assume condition ( 17.22 ) and normalize functions fjk on all oblique edges by a 
condition 

Then, matrices *4i OC ai in all vertices of triangular lattice have the form 



a i 
Aocal = I d 1 g 
h k 



(17.24) 



Proof follows immediately from considering relation (17.23) in the point Pi, taking into account 
the fact that function fjk equals zero in that point if 2 or 4 is contained among the numbers j and 
k, while otherwise equals unity. Note that we do not use conditions (17.21) here. 



It follows from the form (17.24) of matrix -4i OC ai that every closed path on triangular lattice 
to which a zero weight ( = product of all matrix elements, associated with that path, of matrices 
■Aiocal) corresponds, goes either entirely through vertical and horizontal edges or is just an oblique 
line (wound around the torus, of course). The determinant I(u, w) ( 13. 1 ) factorizes in a product of 
a factor corresponding to paths through the square lattice and depending on a, c, h, k and, certainly, 
on u and w, and a factor corresponding to oblique paths and depending only on u and w. 

It follows from the table on p. [55| together with (17.22) that functions u and w on Tq = CP\ 

( Qi C 

have each one zero and one pole. The matrix I , , 
relations ( 14. 6| ), which in our case can be written as 



is obtained from a relation like the first of 



a c 
h k 



fl2 
flS 



whence 



a — u c 
h k — w 



Ufl2 
wf 2 3 



= 0. 



(17.25) 



We will simplify our task still more by considering the "usual" 6-vertex model, assuming con- 
ditions 

a = k. c = h, c 2 — a 2 = 1 



(see Section O, the text between Lemmas 13.1 and 13.2, on the connection between Boltzmann 



weights and matrix elements). Then (17.25) turns into 

aw + 1 



w — a 



(17.26) 



Let us assume that uq and wq in ( 17. lq ) are real and positive. For concreteness, let uo and wq 
be not very different from unity, and 

-1 < a < 1. 



The integration domain A in ( 17.16| ) for thi s cas e (a lune between two circumferences) is depicted, 
in the complex plane of variable w, in Fig. 3.11. 
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Figure 3.11. Integration domain A in the double integral in formula (17.16) in case of a homoge- 
neous model on square lattice, depicted in the complex plane of variable w 



Transform the double integral in ( 17.16] ) into a single integral. We have 



d In u A d In w = 



A 





u 


// dl ° 






u 



Here we used the fact that In 



In 

u 
u 



dlniu = 



In 



-ABC 



A d In w = 



din? 



(17.27) 

is a single-valued function, moreover equalling zero in the 



arc CD A] dA is the boundary of A, i.e. DC ABC. Divide the int egral (J1T.27 ) by as required 
in formula (17.16), and transform it further, using, in particular, (17.26j): 



1 

47ri 

—ABC 

where <po found from conditions 



hi 



4>o 



d In w = 



In 



awoe 



irf> 



1 



-00 



< 00 < 7T, 



awoe 



ua(woe 1 ^ — a) 



uo{woe 1 ^ — a) 



= 1. 



(17.28) 



The integral ( 17.28 ) is the essential, i.e. depending on a, part of the limit in the LHS of ( 17.16 ). 
Of course, it coincides, to within an additive constant (to be exact, a function depending only on 
uq and wq), with the specific free energy of the homogeneous "free fermionic" 6- vertex model on 
the square lattice, see e.g. references j76|, 



58 



18. A "local" approach to orthogonal matrices 



In this section we will show what symmetry conditions must be imposed on the "local" algebro- 
geometrical objects from Section [L4| in order to obtain the description of evolution of orthogonal 



matrices. In contrast to direct calculations of Section 11, here we simply present constructions 
made in the spirit of works |64|, |(3^, |(3(| |6^], and prove their necessary properties. 

Let, as in Section [lj, there be given an algebraic curve T and six divisors in it consisting each of 
one point, namely Pi, . . . , Pg. Let, besides, the curve r , like the "global" curve T from Section |ll| , 
possess an involution under which P±,...,Pq are transformed as follows: Pi <-> P§, P2 <-> P4, 
P3 <-> P5. Consider meromorphic differentials ip on Tq possessing the following properties: a) tp 
has poles (of not higher order than 1) only in points Pi, . . . ,Pg and b) tp transforms into — ip under 
the action of involution I. It follows from those properties and from the fact that the dimension 
of the linear space of holomorphic differentials on Tq equals go (the genus of the curve) that the 
dimension of linear space of differentials ip equals go + 3. 

Fix now a concrete differential ip with non-zero residues in all points Pi , . . . , Pq, and consider 
a divisor T> on To such that V + V 1 is exactly the zero divisor of ip (here V 1 is the image of V 
under involution /). Thus, 



V + V 1 



P 



P fi 



(18.1) 



where P C an is any canonical divisor on To- The relation (18.1) is a "local" analog of relation (11.14). 
Divisor T> consists of go + 2 points, while the dimension of linear space of meromorphic functions / 
satisfying condition 

(f)+V>0, (18.2) 

generically, equals 3. 



Define the scalar product (/, g) of two functions satisfying ( 18.2 ) as the sum of residues of 
differential fg 1 ^ in points Pi, Pi, P5 (g 1 is, of course, the image of g under evolution /): 



(/,<?>= E Res /ffV- 

Pi,-p 2 ,p 5 



(18.3) 



We will study the properties of this scalar product step by step. First let us verify that it is 
symmetric, i.e. that 

(f,9) = {9,f)- (18.4) 

Apply involution / to all objects entering in RHS of ( |l8.3|) , which will cause no change to that 
RHS as a whole. We get: 



if, 9)= E Res(-/V)= E Res /VA- 

P%,P4,,P 3 Pl,P2,P 5 



(18.5) 



The right-hand equality in ( 18.5 ) follows from the fact that the sum of differentia l f I g ip residues, 
taken over all it s pol es Pi, . . . ,Pq, equals zero. Further, the rightmost side of ( 18.5 ) obviously 
equals (g, f), so ( 18.4 ) is proved. 

Consider now scalar products of functions fjk satisfying conditions 



f jk +V-P j -P k >0. 



(18.6) 



In other words, function fjk has a pole divisor T> and must have zeros in points Pj and We 
will prove the equalities 

(18.7) 
(18.8) 
(18.9) 

Let us prove, for example, that (/26;/36) = 0- We have f^ 6 — /51, so meromorphic differential 



(/26, /36) 


= (/56, /46) 


= 0, 


(/l2, /23) 


= (/26, /52) 


-0, 


(/l3, /34) 


= (/36, /23) 


= 0. 



j 26 fio^ has no poles at all in points Pi, P2, P5, whence the scalar product (18.3) indeed equals 0. 
The remaining equalities (18.7 18. 9| ) are proved similarly. 
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If now we manage to normalize all functions fjj. entering in equalities (18.7- 187^) so that 

(Wi*) = l, (18.10) 



and turn to relations (14.5, 14. E) defining "local" matrices 



for I = 1,2,3, we will see 



that those matrices are orthogonal ((18.7), (18.8) and (18.9) are responsible, together with (18.10), 
for orthogonality of local matrices with I = 1, 2 and 3 respectively). At the same time, (18.10) 
will show the non-degeneracy of scalar product (18.3). Thus consider, e.g., a scalar square of 
function /12. We have f[ 2 = /64, so there remains only a residue in point P5 in formula (18.3). 
Generically, /i2(-Ps) f&iiPs) ^ 0, so this residue does not vanish, and we can divide /12 by its 
square root. 

Similarly, the scalar squares of other functions entering in ( |18.7| - p"8~S| ) are non-zero, hence those 
functions can be normalized according to (18.10). 

ai k 



As a result, we have proved so far the orthogonality of three matrices 



Q di 



defined by 



relations ( 14. 5j [14.6 ) and corresponding to vertice s of some chosen triangle in kagome lattice, de- 
picted in Fig.|3.9|(p. p|). To all (nine) ed ges in Fig. |3.9| meromorphic functions corresponded whose 
divisors satisfied equality ( 18.2 ). We ha ve introduced a scalar product in the three- dimen sional 
space of such functions by formula ( |l8.3[) , and normalized the functio ns ac cording to ( 18.10 ). 

We can deal similarly with any other triangle of the form of Fig. 3J3 . If such a triangle is £ 
horizontal lattice periods and 77 vertical periods far from the one that we have considered, then 
divisor T> must be replaced by 



V 



= V + Z(P 1 +P 2 -P i -P 6 )+r 1 (P 3 + P 4 -P 5 - P 2 ). 



Remarkably, the scalar product in the space of meromorphic functions / such that 

(/) + > 



may (and must) be introduced again by formula (18.3) with the same differential tp as for functions 
satisfying ( 18.2 ). The point is that "superfluous" poles of function / are exactly compensated by 
zeros of function g 1 and vice versa, so that differential fg 1 ^ can, as before, have only first order 
poles in points P\, . . . , P§. 



The analogs of orthogonality relations ( 18.7 - 18.9 ) for all lattice vertices are proved by as before. 
In the same way it is proved that we can normalize functions on all edges by setting their scalar 
square equal to unity. The "outer" edges of any diagram of type as in Fig. |3.9| enter also in one 
of neighboring diagrams. It is important that the normalizing of a function on such edge does not 
depend on a diagram to which we consider the edge belongs. 

We will formulate the results of this section as a following theorem. 

Theorem 18.1. If a curve Tq in the situation of Section possesses an involut ion I mapping the 
points Pi,P2,P3 in Pq,P4,Pc, respectively, while a divisor T> satisfie s condition ( 18.1 ), then under 
a proper gauge all square matrices in LHS's of equalities ( 14.4 - 141%) are orthogonal. 



Proof. It remains to give two simple clari ficatio ns. First, the gauge mentioned in the theorem 
is fixed by normalization conditions of type ( 18. 1C ). Second, we considered in this section only a 
situation at one moment of discre te ti me. H owever, a simple calculation based on divisor evolution 
described between Lemmas 



14.1 



and 



14.2 



shows that our considerations are valid for any time 



("superfluous" poles of each function in ( 18.3 ) are compensated by zeros of another function as 
well as before). 



19. Reduction to Ising model 

It is known that Onsager's "star-triangle" transformation (|33j, see also manual |76|), which can 
be graphically represented as 

X - A, (19.1) 



GO 



converts the statistical mechanical Ising model on a plane hexagonal lattice to the Ising model on 
a triangular lattice. Imagine now the obtained triangular lattice as made up of triangles of the 
form \/ . If we apply the triangle-star 

V+Y, (19.2) 

transformation to those triangles, we will come to hexagonal lattice again. If Ising model is 
inhomogeneous, i.e. a coupling along a given edge depends on the edge, then transformations 



( 19.1 ) and ( 19.2 ) applied alternately lead to some evolution of those coefficients on alternating 
hexagonal and triangular lattices. A hypothesis looks natural that this evolution is completely 
integrable. 

In this, not very large, section we will show that the mentioned evolution of Ising model coupling 
coefficients is isomorphic to the evolution of orthogonal matrices in vertices of triangular lattice. 
The latter evolution was considered from a "local" algebro-geometrical viewpoint in Section [l8| 

The direct way to Ising model turns out to go via considering the orthogonal matrices with 
determinant —1. As we will see, this is connected with the fact that such matrices are symmetric 
(recall that they are of size 2x2). 

So, consider a relation 



a i 


h 




h 


-CLi 


•) 


















(i 




K 




b' 3 


-a' 3 










A3 




&3 


-a 3 




b't 


K 


-a[ 









(19.3) 



where for all j 



a ) + q = i, (a> J f + (b>f = i. 



The relation ( 19. 3| ) corresponds to description of matrix evolutio n give n in Section |TJ: at every 



step a 3 x 3 matrix obtained as a product of the form as in RHS of ( 19.3 ) is factorized in a product 
of the form as in the LHS of that formula. 

Introduce Ising model coupling coefficients Ki, K 2 , K 3 , L\, L 2 , L 3 by formulae 

e^p(±2K 1 ) = 61 T iai, exp(±2 J ftT 2 ) = ^ 2 T ^ , 

exp(±2i4r 3 ) = 63 T ia 3 , (19.4) 



exp(±2i!) = — ^ — , exp(±2L 2 ) = b' 2 =p iaL 

exp(±2L 3 ) = !^i±l)_ (19 5) 

a 3 

The following lemma is principal in this section. 

Lemma 19.1. If ( jl9.4 ) holds, the coefficients K\, . . . , L3 satisfy the following star-triangle rela- 
tions (we reproduce formulae (6.4.8a-d) from the book ydj ; coefficients Kj belong to a triangle, 



Lj — to a star, R is some numeric factor): 

2cosh{L 1 +L 2 +L 3 ) = Rexp{ K i+K 2 + K 3 ), (19.6) 

2cosh(-L 1 +L 2 +L 3 ) = Rexp(K 1 -K 2 -K 3 ), (19.7) 

2cosh(ii-L2 + L 3 ) = Rexp(-K 1 + K 2 -K 3 ), (19.8) 

2 cosh(ii + L 2 - L 3 ) = i?exp(-Xi- K 2 + K 3 ). (19.9) 
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Figure 3.12. Star-triangle transformation is shown by thick lines, while the corresponding matrix 
transformation — by thin lines 




Figure 3.13. Triangle-star transformation is shown by thick lines, while the corresponding matrix 
transformation — by thin lines 



Lemma 19.1 is illustrated by Fig. 3.12 



Proof of Lemma 19.1 will be done using the following parameterization of relations ( 19. 3| ) with 
variables k, A, /i, where k is the modulus of all elliptic functions: 



snA, a2 = fcsn (A + /i), <23=sn/i, 



(19.10) 



6i=cnA, fe 2 = dn(A + /i), &3=cn/x, (19.11) 
a' 1 = ksnX, a 2 =sn(A + /i), a 3 = fcsn//, (19.12) 
&i=dnA, fr 2 = cn(A + ^), b' 3 = dnfi. (19.13) 
The fact that ( [l 9 . 1 C| — p~9 . 1 3| ) is really a parameterization of ( 19.3 ) is verified directly using elemen- 



tary properties of elliptic functions. Then, from the same properties it follows that formulae for 
K 2 and L 2 from ( |l9.4[ ) and (|l9.5| ) may be rewritten as follows: 



where 



exp(±2X2) = cn k =F *sn k, exp(±2L2 



A + k + /J, = il . 



z(dn k ± 1) 
k sn k ' 



(19.14) 



(19.15) 



Here we have denoted the half-period of elliptic functions by following R. Baxter's book |76|], 
although it is usually denoted K' . Our formulas (19.4) and (19.5) coincide now, with regard to 
(19.14), with parameterization (7.8.5) from |7(| for relations (19.6 - 19j| ), while ( 19.15Q coincides 
with formula (7.13.4) from ]7q| if our X,k,/i equal Baxter's iui,iii2,iu3. The lemma is proved. 



We have described the transformation (19.1) by means of orthogonal matrices, and illustrated 
this with Figure |3.12 Similarly, transformation (19.2) together with the transformation of orthog- 
onal matrices is illustrated by Figure 3.13. The latter differs from Fig. 3.12 in the fact that now 



G2 



V 



to a kagome lattice triangle 



corresponds an Ising triangle, while to a triangle 



corresponds a star. One can reduce this situation to the previous one by transposing the products 
of orthogonal matrices (recall that matrices themselves are symmetric), corresponding to LHS and 
RHS of Fig. 3.13| and then changing their places. Then it is clear that relations of the form (19.4. 
19.5) between matrix elements and coefficients Kj and Lj for triangle and star respectively work 
in the case of Fig. 3.13 as well. 
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